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ABSTRACT 

The magnetar IE 1547.0-5408 exhibited outbursts in October 2008 and January 2009. In this 
paper we present in great detail the evolution of the temporal and spectral characteristics of 
the persistent total and pulsed emission of IE 1547.0-5408 between ~1 and 300 keV starting in 
October 3, 2008, and ending in January 2011. We analyzed data collected with the Rossi X-ray 
Timing Explorer, the International Gamma-Ray Astrophysics Laboratory and the Swift satellite. 
We report the evolution of the pulse frequency, and the measurement at the time of the onset of 
the January 2009 outburst of an insignificant jump in frequency, but a major frequency derivative 
jump Ai> of +(1.30 ± 0.14) ■ lO" 11 Hz/s (Az>/z> of -0.69 ± 0.07). Before this v glitch, a single 
broad pulse is detected, mainly for energies below ~ 10 keV. Surprisingly, ~ 11 days after the 
glitch a new transient high-energy (up to ~ 150 keV) pulse appears with a Gaussian shape and 
width 0.23, shifted in phase by ~ 0.31 compared to the low-energy pulse, which smoothly fades 
to undetectable levels in ~ 350 days. We report the evolution of the pulsed emission spectra. 
For energies 2.5-10 keV all pulsed spectra are very soft with photon indices T between -4.6 and 
-3.9. For ~ 10-150 keV, after the v glitch, we report hard non-thermal pulsed spectra, similar 
to what has been reported for the persistent pulsed emission of some Anomalous X-ray Pulsars. 
This pulsed hard X-ray emission reached maximal luminosity 70±30 days after the glitch epoch, 
followed by a gradual decrease by more than a factor 10 over ^300 days. These characteristics 
differ from those of the total emission. Both, the total soft X-ray (1-10 keV) and hard X-ray (10- 
150 keV) fluxes were maximal already two days after the 2009 January outburst, and decayed by 
a factor of >3 over ~400 days. The total spectra can be described with a black-body (kT values 
varying in the range 0.57-0.74 keV) plus a single power-law model. The photon index exhibited 
a hardening (~-1.4 to ~-0.9) with time, correlated with a decrease in flux in the 20-300 keV 
band. We discuss these findings in the framework of the magnetar model. 

Subject headings: pulsars: individual (IE 1547.0-5408; SGR J1550-5418; 1RXS J1708-4009; 
4U 0142 + 61; IE 1841-045), X-rays: stars 
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1. Introduction 

Over the last decade, observational evidence 
has mounted that Anomalous X-ray pulsars 
(AXPs) and Soft Gamma-ray Repeaters (SGRs) 
belong to the same class of objects. Most publi- 
cations discuss these objects in the frame work of 
strongl y magnetized neutron stars, dubbed mag- 



netars ( Duncan fc Thompson! 1992 ; Thompson fc Duncan 
1995). Alternative interpretations exist including 
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'normal' m agnetized neutron stars with fall-bac k 
disks (e.g. IChatteriee et all l2000) : lAlpar.1 [iool. 



white dwarfs (e.g^ Malhcir o et al 



quark stars (e.g. lXu 1120071; lOrsaria et alj|201lh . In 



2011), or even 



this paper we adopt the magnetar interpretation. 
Magnetars are rotating neutron stars with surface 
magnetic field strengths of 10 14 — 10 15 G, well 
above the critical field strength of 4.413 x 10 13 
G at which the cyclotron energy of an electron 
reaches the electron rest mass energy. They have 
long rotation periods P (in the range 2-12 s) and 
large period derivatives P (~ 10~ n s s _1 ). 

Their soft X-ray (<10 keV) luminosities of ~ 
10 erg s _1 exceed the available rotational en- 



ergy l osses (see reviews by IWoods fc Thompson 
2006t iKaspil l2007t iMereghettil l2008h . In addi- 
tion to SGR giant flares, SGRs and AXPs exhibit 
periods of high activity wherein radiative out- 
bursts are accompanied by numerous short bursts 
of typical durations a few hundred milliseconds. 
In the context of the magnetar twisted magne- 
tosph ere model ([Thompson. Lvutikov fc Kulkarni 
20021 ). the energy release is due to magnetic field 
re-arrangement which may be triggered by crustal 
deformation causing a glitch in their rotational 
timing behaviour. 

A distinct feature is the high non-thermal 
luminosity of the persistent emiss ion of AXPs 



above 20 keV, first esthablished by iKuiper et al 



(|2004l ) with the discovery of pulsed hard X-ray 
emission up to 150 keV from AXP IE 1841- 
045 in Supernova remnant (SNR) Kes 73. The 
X-ray luminosity spectrum exhibits two peaks, 
one near 1 keV and the other above 100 keV 



(see e.g. Kuioer et al. 


20061: G6tz et al. 20061: 


den Hartog et all 2008a 


b; Enoto et al. l2010bh. 



The soft X-ray part of the spectrum below 10 keV 
and the hard X-ray part above 20 keV, both ex- 
hibit luminosities exceeding the spin-down power 
1-3 orders of magnitude. The spectra of the 
first component can be empirically described by 
a thermal black body (BB) plus a soft power- 
law (index 2-4) model or two BB models, and 
more phy sically with resonant cyclotron sca t tering 



models (Thompson. Lvutikov & Kulkarni 200 



Lvutikov fe Gavrhll 2006 ^Fernandez fc Thompson! 



2007t iGiiver et alJ|2007t l2008t iRea et all 120081 
Nobili et al.ll2008t fzane et al.ll2009l l201ll) . 

On the origin of the hard non-thermal per- 
sistent component above ~ 10 keV, however, 



no con census has been reache d. It has been 
shown ( den Hartog et al. 2008al lb) that this com- 
ponent is persistent and within the statistical 
errors of about 20% stable over years, possi- 
bly as long as a decade, in both total and 
pulsed fluxes and spectral shapes, as well as in 
pulse phase/shape. Most recent attempts to ex- 
plain these fi ndings in the magnetar i n terpre - 



tation are by iBeloborodov fc Thompson! (1200 
Beloborodovf (|2009|): tearing fc Harding) ([200 



2008); lPavan et al 



7); 
07. 



20091 ). The non-thermal com- 



ponent above 10 keV has also been discussed in the 



context of the fall-back disk model (Triimpcr et al. 
2010al lbh. In his latest work, beloborodovl (12009. 



20111 ) discusses how a "starquake" can cause con- 
vective motions in the crust which twist the mag- 
netic field anchored to the surface, after which it 
gradually untwists, dissipating magnetic energy 
and producing radiation. A current carrying bun- 
dle of closed field lines ("j-bundle") is created 
above the magnetic dipole axis. Part of the ther- 
mal X-rays emitted by the neutron star are up- 
scattered by the inner relativistic outflow in the 
quasi-steady j-bundle, producing a beam of hard 
X-rays. 

In this context, it is therefore of great inter- 
est to study the evolution of the observational 
characteristics of a magnetar before and after a 
glitch followed by a radiative outburst; to study 
the spectral variations in the persistent steady 
and pulsed emissions in the soft and hard X- 
ray bands, as well as variations in pulse profile. 
Sofar the decay from a radiative outburst could 
only be studied in detail for energies below ~ 10 



keV, e.g. for 1E2 259+586 (Woods et all 1200 



SGR 0501+4516 (R ea et al 



2009; Gogiis et al. 



2010al) . SGR J1 833-0832 dGogiis et all 120101 
SGR 1900+14 dGogiis et all 120111). SGR1833- 
0832 (lEsposito et alJl201lh a nd CXOU J164710.2- 



455216 (|Woods et all l201lh . In this work we 



present the first study in which both the tem- 
poral and spectral characteristics at both soft and 
hard X-ray energies before and during a radiative 
outburst of a magnetar, IE 1547.0-5408, have been 
derived. 

The X-ray source IE 1 5 47.0-5 408 was discov- 
ered by lLamb fc Marked (jl98ll ) in an Einstein 
HRI observation Based on its X-ray spec trum 
and variability, iGelfand fc Gaenslerl (j2007) dis- 
cussed this source as candidate magnetar in a can- 
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didate supernova remnant (G327.24-0.13). Sub- 
sequently, it s period was discov ered in the ra- 
dio band by ICamilo et all (|2007h . making it the 
first magnetar of which the pulsation was not dis- 
covered at X-ray energies, and the second known 
radio-emittin g magnetar after t he transient XTE 



J1810-197 (jCamilo et all l2006h . With the re 



ported period of 2.069 s it appeared to be the 
fastest spinning magnetar known, with a derived 
surface magnetic dipole field strength of 2.2 x 10 14 
G, a characteristic age of 1.4 kyr, and a spin-down 
luminosity of 1.0 x 10 35 ergs s _1 , assuming a dis- 
tance of ^9 kpc (estimated from the dispersion 
measure). X-ray puls ations were fir s t det ected 
with XMM- Newton by lHalpern et~ai1 (j2008l) in a 
period of enhanced activity in 2007. These authors 
concluded that the source was recovering from an 
X-ray outburst in between 2006 and 2007. 

On Oct. 3, 2008 (MJD 54742) IE 1547.0-5408 
started a perio d of strong burstin g activity de- 



tecte d by Swift (jKrimm et a l. 2008a.bl; llsrael et al 
2010 ), and by the Gamma-ray Burst Monitor 



(GBM) on the Fermi Gamma-ray Space Telescope 
|Kaneko et al.ll2010l) . By the latter authors the 
source was dubbed SGR J1550-5418, motivated 
by the hundreds of typical soft-gamma-repeater 
bursts. The bursting activity was accompanied by 
a radiative outburst which was also d etected with 
RXTE and Chandra dNg et alj|201lh. Ilsrael et all 
(|2010h used Swift and lNg et al] (|201lh RXTE and 
Chandra data to study the evolution of the pulse 
profile shape and the phase-averaged spectrum for 
energies below 10 keV during the outbursting ac- 
tivity. 

The Swift Burst Alert Telescope (BAT) re- 
ported renewed extreme bursting activity starting 



on J an. 22, 2009 (MJD 54853) (Gronwall ct al 



2009), also detected by I NTEGRAL (jSavchenko et al 



20091: iMereghetti et alj|2009al) . the Fermi GBM 



also part of our 600 ks INTEGRAL open-time 
(PI den Hartog) ToO observations, a power-law 
index of ~ 1.5 was found for energies from 20 
up to 150 keV from 2 to 7 days afte r the Jan. 
2009 outburst (|den Hartog et al.ll2009l) . This was 
confirmed in Suza ku observations tak en 7 days 
after t he outburst (lEnoto et al. 2010al ). Further- 



Kuiper et al.l (|2009al ) announced the de- 



more, 

tection with INTEGRAL of pulsed emission for 
energies up to 150 keV. 



Enoto et al.l (|2010al ) presented the Suzaku re- 



sults on AXP IE 1547.0-5408 for their observa- 
tions 7 days after the onset of the bursting activ- 
ity in 2009. Pulsations were detected for energies 
up to 70 keV, and the total, ti me-averag e d 0.7 - 
114 keV spectrum was discussed. Ng et al.l ( 201ll ) 
compared the 2008 and 2009 outbursts, analyzing 
RXTE and Chandra observations taken during ~ 
20 days after the onset of the bursting activities in 
2008 and 2009, comparing pulse profiles, pulsed- 
emission spect ra below 10 keV, a nd spi n evolution. 
Furthermore, iBernardini et al.l (|201ll ) monitored 
the 2009 outburst with Chandra, XMM-Newton 
and INTEGRAL over a period of two weeks after 
the onset, as we ll as with Swift over a 1.5 year 
interval. Finally, Scholz fc Kaspi I (|201ll ) focussed 
on the 2009 outburst using Swift XRT observa- 
tions for energies below 10 keV addressing the per- 
sistent radiative evolution and a statistical study 
of the burst properties. 

In this work we exploited the extensive data 
base of observations of IE 1547.0-5408 with RXTE, 
INTEGRAL and Swift covering the 27 months 
from the onset of the bursting activities in Oc- 
tober 2008 till January 2011. We studied in 
detail the variable (total and pulsed) emission, 
but did not address the burst properties (see 
e.g. ISavchenko et alJlioiot IMereghetti et al.ll2009t 
Kan eko et al.ll2010l ). The evolution of the tim- 



2009i1. RHESSI (|BeUm et al.l l2009h. a nd the Wide- 



band All-sky Monitor aboard Suzaku (|Terada et al 
20091) . Follow-up observations with Swift and 



( Connaughton fc B rigs _ ^onus-Wind (jGolenetskii ein alparameters. pulse profiles, fluxes, spectra have 

been derived over the broad energy range ~ 1 - 
300 keV. Particularly interesting is the discovery 
of a transient non-thermal hard X-ray component, 
which appeared shortly after the onset of Jan. 22, 
2009 outburst, triggered by a v timing glitch, as a 
distinct new pulse in the pulse profile. Finally, our 
findings are summarized in Sect. [5] and are com- 
pared with model predictions within the magnetar 
framework. 



XMM-Newton revealed dust scattering X-ray 
rings centered on IE 1547.0-5408, from which 
Tiengo et al. (2010) deduced a source distance 
of ^3.9 kpc. Based on a 100 ks public Target- 
of-Opportunity (ToO ) INTEGRAL observation 
Baldovin et al. ( 20091) reported the discovery of a 
hard power-law tail in the spectrum. Exploiting 
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2. Instruments and observations 
2.1. Rossi X-ray Timing Explorer 

In this study extensive use is made of data from 
monitoring observations of AXPs with the two 
non-imaging X-ray instruments aboard RXTE, 
the Proportional Counter Array (PCA; 2-60 keV) 
and the High Energy X-ray Timing Experiment 
(HEXTE; 15-250 keV). 

2.1.1. RXTE PCA 

The PCA (jjahoda et all Il996l) consists of 



Table 1: List of RXTE observations of IE 1547.0- 
5408 used in this study. 



five collimated xenon proportional counter units 
(PCUs) with a total effective area of ~ 6500 
cm 2 over a ~ 1° (FWHM) field of view. Each 
PCU has a front Propane anti-coincidence layer 
and three Xenon layers which provide the basic 
scientific data, and is sensitive to photons with 
energies in the range 2-60 keV. The energy res- 
olution is about 18% at 6 keV. All PCA data 
used in this study have been collected from ob- 
servations in GoodXenon mode allowing high-time 
resolution (0.9/xs) analyses in 256 spectral bins. 
Since the launch of RXTE on Dec. 30, 1995 the 
PCA has experienced high voltage breakdowns for 
all constituting PCUs at irregular times. To avoid 
further breakdowns, during its already 14.5 year 
mission not all PCUs are simultaneously operat- 
ing. The most stable PCU is PCU-2, which is on 
for almost all of the time. On average one (50%) 
or two (40%) PCUs is/are operational during a 
typical observation. 

2.1.2. RXTE HEXTE 



The HEXTE instrument (Rothschild et al 



19981 ) consists of two independent detector clus- 
ters A& B, each containing four Na(Tl)/ CsI(Na) 
scintillation detectors. The HEXTE detectors are 
mechanically collimated to a - 1° (FWHM) field 
of view and cover the 15-250 keV energy range 
with an energy resolution of ~ 15% at 60 keV. 
The collecting area is 1400 cm 2 taking into ac- 
count the loss of the spectral capabilities of one 
of the detectors. The maximum time resolution 
of the tagged events is 7.6/xs. In its default op- 
eration mode the field of view of each cluster is 
switched on and off source to provide instanta- 
neous background measurements. However, also 
HEXTE suffers from aging, and since July 13, 



Observation 


Begin/End Date 


Exp.W 


identifier 


(yyyy/ 


mm/dd) 


(ks) 


20060 


1997-05-15 


1997-05-15 


0.376W 


93017 


2008-10-03 


2009-01-24 


162.484 


94017 


2009-01-25 


2009-06-22 


174.676 


94427 


2009-06-30 


2009-12-26 


104.677 


95312 


2010-01-04 


2010-12-25 


236.952 


93017 - 95312 


2008-10-03 


2010-12-25 


678.789 



l 2 l PCU-2 exposure after screening 



2006 HEXTE cluster-A operates in staring mode 
at an on-source position, while on March 29, 2010 
cluster-B was commanded to stare at an off-source 
position. 

Due to the co-alignment of HEXTE and the 
PCA, they simultaneously observe celestial tar- 
gets. Table Q] lists the RXTE observations used in 
this study with in the fourth column the screened 
exposure of PCU-2 (see Sect[3|). A typical obser- 
vation consists of several sub-observations spaced 
more or less uniformly (4 — 9 days apart) between 
the start and end date of the observation cycle. 
However, for a ~ 2 — 3 weeks time period directly 
after the two outbursts intensive timing has been 
performed. 

2.2. INTEGRAL 



The INTEGRAL spacecraft (jWinkler et al 
2003), launched 17 October 2002, carries two 



main 7-ray ins truments: a high-ang ular-resolution 



imager IBIS (jUbertini et al.l 120031 ) and a high 



energ y-resolution spectrometer SPI (jVedrenne et al 
20031 ). The pay load is further supported by two 
monitor instruments providing complementary ob- 
servations in the X-ray and optical energy bands, 
the Joint European Monitor for X-rays (JEM- 
X; iLund et al. 2003) and the Optical Monitoring 
Camera (OMC; 500 - 600 nm (Johnson V-filter)). 
All the high-energy instruments make use of coded 
aperture masks enabling image reconstruction in 
the hard X-ray/soft 7-ray band. 

In our study, guided by sensitivity consid- 
erations, we only used data recorded by the 
INTEGRAL Soft Gamma-Ray Imager ISGRI 
(jLebrun et al.l 120031) , ;he upper detector system 



4 



Table 2: List of INTEGRAL observations of 
IE 1547.0-5408 used in this work, sorted on INTE- 
GRAL orbital revolutions. In the last two columns 
are given the exposures and the time segment iden- 
tifiers. 



Revs. Begin/End Date Begin/End 

(yyyy/mm/dd) m.id 



Exp. I 1 ! 
(kg) 



731 

767 
768-772 
782-791 
840-850 
899-910 
911-912 



2008- 10-08 

2009- 01-24 
2009-01-28 
2009-03-09 

2009- 08-30 

2010- 02-23 
2010-03-30 



/ 2008-10-10 
/ 2009-01-25 
/ 2009-02-08 
/ 2009-04-08 
/ 2009-10-01 
/ 2010-03-29 
/ 2010-04-03 



54747.9 - 54749.2 
54855.6 - 54856.9 
54859.6 - 54870.7 
54899.6 - 54929.0 
55073.6 - 55105.5 
55250.3 - 55284.3 
55285.9 - 55289.0 



98.4 
92.6 
575.6 
191.9 
407.8 
524.0 
193.7 



^'Effective on-axis exposure 



of IBIS, sensitive to photons with energies in the 
range ~20 keV - 1 MeV, and JEM-X operating in 
the 3-35 keV X-ray band. 

With an angular resolution of about 12' and a 
source location accuracy of better than 1' (for a 
> lOcr source) ISGRI is able to locate and sepa- 
rate high-energy sources in crowded helds within 
its 19° x 19° field of view (50% partially coded) 
with an unprecedented sensitivity (~ 960 cm at 
50 keV). Its energy resolution of about 7% at 100 
keV is amply sufficient to determine the (contin- 
uum) spectral properties of hard X-ray sources in 
the ~ 20 - 300 keV energy band. 

The timing accuracy of the ISGRI time stamps 
recorded on board is about 61/^s. The time 
alignment between INTEGRAL and RXTE is 
better than ~ 50/^s, verified using data from 
simultaneous RXTE and INTEGRAL observa- 
tions of the accretio n-powered milliseco nd pulsars 
IGR J00291 +5934 (iFalanea et alJ l2005l) and IGR 
J17511-3057 (|Falanga et alj|201lh . 

JEM-X consists of two identical telescopes each 
having a field of view of 7? 5 (diameter) at half 
response and able to pinpoint a 15cr source with 
a 90% location accuracy of about 1'. Its energy 
resolution and timing accuracy (3cr) are 1.3 keV 
at 10 keV and 122. 1/xs, respectively. 

In its default operation mode INTEGRAL ob- 
serves the sky in a dither pattern with 2° steps, 
which could be rectangular e.g. a 5 x 5 dither 
pattern with 25 grid points, or hexagonal with 7 
grid points (target in the middle). Typical inte- 



id. 



gration times for each grid point (pointing/sub- 
observation) are in the range 1800 - 3600 seconds. 
This strategy drives the structure of the INTE- 
GRAL data archive which is organised in so-called 
science windows (Sew) per INTEGRAL orbital 
revolution (lasting for about 3 days) containing 
the data from all instruments for a given point- 
ing. Most of the INTEGRAL data reduction in 
this study was performed with the Offline Scien- 
tific Analysis (OSA) version 7.0 distributed by the 
INTEGRAL Science D ata Centre (ISDC; see e.g. 
Courvoisier et al.1l2003l ). 

Table lists the INTEGRAL orbital rev- 
olution (Rev.) identifiers with corresponding 
start/end dates of the observations used in the 
imaging/spectral analyses and/or timing analyses 
of IE 1547.0-5408. Revolutions 731 (Oct. 2008; 
ToO), 767 (Jan. 2009; public ToO), 768-772 (Jan- 
Feb., 2009; ToO) and 911-912 had IE 1547.0-5408 
as on-axis prime target, while during revolutions 
782-791, 840-845/850 and 899-910 (all part of an 
INTEGRAL key program to observe deeply the 
Z=±25 Galactic plane region) IE 1547.0-5408 was 
in the field of view at moderate off-axis angles. 

2.3. Swift 



The Swift satellite (jGehrels et al.l 120041 ) was 
launched on November 20, 2004, and regular ob- 
servations began on April 5, 2005. The main goal 
of the mission is the study of gamma-ray bursts 
and their afterglows. Swift carries three co-aligned 
instruments: The wide-field coded aperture mask 
Burst Alert Telescope (BAT; 15-150 keV), the nar- 
row field (23f 6 x 23.'6) grazing incidence Wolter 1 
X-Ray Telescope (XRT; 0.2-10 keV) and the Ul- 
traviolet/Optical Telescope (UVOT). 

During periods of non-burst operations obser- 
vations of other high-energy sources are scheduled. 
In this work we us ed data of IE 1547.0 -5408 gath- 
ered by the XRT (jBurrows et alj|2005l) . 

Regular monitoring of IE 1547.0-5408 com- 
menced on June 22, 2007 following the detection 



of ra dio pulsations on June 8, 2007 (jCamilo et al 
2007). The XRT has several operation modes 



of which we used only the Windowed Timing 
(WT) and Photon-counting (PC) modes. In PC 
mode full imaging (600 x 600 pixels; pixel scale 
2'.'36/pixel) and spectroscopic resolution are re- 
tained, but its time resolution is only 2.5073 s, 
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insufficiently to detect the pulsations of IE 1547.0- 
5408. In WT mode 10 rows are compressed in one, 
and only the central 200 columns are read out. 
Therefore only a ~ 8' central strip of the field of 
view is covered. The time resolution in WT mode, 
however, is 1.7675 ms amply sufficient for pulse 
timing studies of IE 1547.0-5408. 

Intensive monitoring data are available for 2- 
3 weeks time periods directly after the onsets of 
the Oct. 2008 and Jan. 2009 outbursts. We 
used the XRT WT data gathered between Oct. 
3, 2008 and Oct. 24, 2008 (Swift observation ids. 
00330353000—00330353016; MJD 54742-54763) 
for a combined RXTE PCA/Swift XRT timing 
study of the pulsed signal at soft X-rays. 

From the onset of the Jan. 2009 outburst on 
Jan. 22, 2009 till Feb. 7, 2009 IE 1547.0-5408 
was observed daily by the XRT both in WT and 
PC mode. This period of dense sampling was fol- 
lowed by monitoring observations taken at a rate 
of roughly 1-4 times per month in either WT or 
PC mode, and this strategy continues up to now. 
XRT WT data gathered between Jan. 22, 2009 
and Feb. 22, 2009 (MJD 54853-54884) have been 
used in the construction of timing models in com- 
bination with RXTE PCA data. 

Data taken in PC mode on Jan. 25 & 29 and 
Feb. 4, 2009 (Swift observation ids. 00341114000, 
00030956034 and 00030956039, respectively) cov- 
ering the time window for which we have deep IN- 
TEGRAL observations (see Table[2] Jan. 24 - Feb. 
8, 2009 i.e. MJD 54855.6-54870.7), have been em- 
ployed in spectral analyses of the total soft X-ray 
emission of IE 1547.0-5408. Furthermore, to study 
the evolution of the total soft X-ray spectrum 
we analyzed Swift XRT observations 00030956048 
(July 8, 2009), 00030956051 (Aug. 19, 2009), 
00090404003 (Apr. 12, 2010), 00090404019 (Sept. 
28, 2010) and 00090404027 (Feb. 25, 2011). 

Finally, XRT observations in WT mode per- 
formed between Jan. 4 and 13, 2009 and after 
Feb. 22, 2009 have been used for verification pur- 
poses of timing models based on merely RXTE 
(monitoring) observations. The latter are sam- 
pled t oo sparsely given t he high level of timing 
noise ([Camilo et al]|2008l ) present in this source. 



3. Timing, total flux and pulsed flux evo- 
lutions below ~ 10 keV 

We exploited the regular monitoring observa- 
tions with the PCA on RXTE over more than 
2 years (see Table []} that provided the required 
high-statistics data. Because the number of active 
PCUs at any time was changing, we treated the 
five PCUs constituting the PCA separately. Good 
time intervals (GTI) have been determined for 
each PCU by including only time periods when the 
PCU in question is on, and during which the point- 
ing direction is within 0?05 from the target, the el- 
evation angle above Earth's horizon is greater than 
5°, a time delay of 30 minutes since the peak of a 
South- Atlantic-Anomaly passage holds, and a low 
background level due to contaminating electrons, 
as measured by PCU-2, is observed. Furthermore, 
periods during which a detector break-down (see 
Sect l2.1.Tj) has occurred, are excluded for further 
analysis. Finally, because we are only interested 
in the persistent non-burst emission properties, we 
rejected small duration (< 1 s) bursts. The lat- 
ter filtering is especially important just after the 
onset of the Oct. 2008 and Jan. 2009 outbursts, 
when copious numbers of bursts are detected. The 
resulting good time intervals have subsequently 
been applied in the screening process to the data 
streams from each of the PCUs (e.g. see Table [1] 
for the resulting screened exposure of PCU-2 per 
observation cycle). 

The TT (Terrestrial Time) arrival times of the 
selected events (for each sub-observation and PCU 
unit) were converted to arrival times at the so- 
lar system barycenter (in Barycentric Dynamical 
Time (TDB) time scale) using 1) the JPL DE200 
solar system ephemeris, 2) the instantaneous 
spacecraft position, and 3) the sub-arcsecond ce- 
lestial position of IE 1547.0-5408. The position 
used is: (a,S) = (15 h 50 m 54fll, - 54°18'23'/7) for 



epoch J2000 (jCamilo et all l2007t ). which corre 



sponds to (l,b)= 
coordinates. 



=(327.23705,-0.13162) in Galactic 



3.1. Total flux evolution at X-rays below 
- 10 keV 

From the screened event datasets and the 
break-down/burst corrected exposure times we 
could derive the total count rate from a circular 
field with a radius of ~ 1° (PCA field of view up to 
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zero collimator response) on IE 1547.0-5408. We 
determined this for the 4-27 PHA channel range 
(~ 2 — 10 keV) including all Xenon detectors 
layers. The count rate includes the instrumen- 
tal and celestial (Galactic diffuse) background, 
emission from IE 1547.0-5408 and its PWN/SNR 
(jVink fc Bambal l2009h and other discrete point 
sources. We found in the archive two RXTE scan- 
ning observations, both performed on May 15, 
1997, with IE 1547.0-5408 in the field of view for 
40 s (screened) at 28.' 1 offset angle and 336 s 
(screened) at 17f3 offset angle. We used the count 
rate from the combination of these 2 scanning ob- 
servations as reference, because the source at that 
time was very likely in quiescence. 

The total flux evolution, from Oct. 3, 2008 up 
to December 25, 2010, of IE 1547.0-5408 above the 
reference level of 9.70 ± 0.16 c/s as measured by 
PCU-2 in the 4-27 PHA range (all detector layers) 
is shown in the upper panel of Fig. [1] 

We cropped this figure in vertical direction in 
order to better visualize the late time evolution 
(the count rate of 44.19(8) measured directly af- 
ter the Jan. 22, 2009 outburst at MJD 54853.911 
is off scale). The onsets of the Oct. 3, 2008 and 
Jan. 22, 2009 outbursts are indicated in this fig- 
ure. There is a monitoring gap of about a month 
which occurred in Sept./Oct. 2009 due to RXTE 
spacecraft anomalies. 

It is clear that since the Oct. 3, 2008 the total 
flux decayed to almost its reference level just be- 
fore the Jan. 22, 2009 outburst. The flux decay of 
the much more intense second outburst, however, 
has still not reached its reference level ~ 2 year 
after the onset. We studied the decay of the to- 
tal flux since the Jan. 22, 2009 outburst in more 
detail by fitting the data, excluding the time in- 
tervals in which the mini-outbursts occurred (see 
Sect l3.1.i|) . with a model composed of a constant, 
an exponential- and a power-law like decay com- 
ponent: 

/ t — t \ a t — t 
Rtot{t)=a + b- 1 + +c-cxp( -) 

V T o / n 

(1) 

We fixed t at MJD= 54853.035, used only times 
with t > 54854 (MJD) and fitted the 6 free pa- 
rameters to obtain the best desciption of the ob- 
served decay. We found the following optimized 
parameters: a = 0.33 ± 0.05 c/s, b = 152.5 ± 4.4 



c/s, c = 3.5 ± 0.1 c/s, r = 0.129 ± 0.004 day, 
n = 207±llday, and finally a = -I.01i0.01 (er- 
rors are ler), for a fit with a reasonable quality of 
xl = 138.58/82 = 1.69. This model is superposed 
in the top panel of Fig. [T] as a solid line, and fol- 
lows the measurements globally pretty well. The 
slowly decaying exponential component has a time 
scale of typically ~ 200 days, whereas the rapidly 
decaying power-law component contributes only 
significantly to the total flux for the first ~ 50 
days since the Jan. 22, 2009 outburst, over which 
its fractional contribution to the total rate is more 
than 10%. The power-law component dominates 
only during a time period of ~ 5 days directly af- 
ter the outburst, then the combination of constant 
plus exponential component takes over. 

2oo | | | | 




54920.58 54920.59 54920.60 54920.61 54920.62 54920.63 
T [MJD;TT] 

Fig. 2.— Count rate of IE 1547.0-5408 in the 3- 
27 PHA range (all detector layers combined) as 
observed by the PCU-2 during observation 94017- 
09-09-00 on Mar. 30, 2009 (MJD 54920). The 
time bin width is 1 s. The time is expressed in 
MJD for the TT time system. Note the enhanced 
emission level of ~ 450 s duration with superposed 
three short bursts. 

3.1.1. Mini- outbursts 

A detailed look at the decaying tail of the Jan. 
22, 2009 outburst shows a large deviation from 
the global decay trend near MJD 54920 (Mar. 
30, 2009). Initially, we believed that we were 
dealing with a processing error/anomaly, but a 
deeper investigation revealed a period of enhanced 
emission of about 450 s with superposed at least 
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Fig. 3.— Count rate of IE 1547.0-5408 for the ob- 
servation 95312-01-02-00 on Jan. 11, 2010 (MJD 
55207) as in Fig. [5] Note the enhanced emis- 
sion level of ~ 130 s duration with superposed 
two bursts separated at 11 s (not resolved at the 
displayed scale). 



three short duration bursts (see Figj2]). This be- 
haviour was seen in both PCU-1 and PCU-2, 
which were both operational during the observa- 
tion in question, thus excluding an instrumental 
(non-celestial) cause of the observed phenomenon. 
The period of enhanced emission was preceded by 
a short duration burst. 

We searched for more mini-outbursts and found 
one on Jan. 11, 2010 (MJD 55207.8845 TT, see 
Fig. [3]). Enhanced emission is seen for a period 
of about 130 s with superposed two bursts with a 
separation of ~ 11 s between the burst maxima. 

An interesting feature of this enhanced emis- 
sion period is that the burst trail is highly pulsed. 
We found a pulsed signal strength of ~ 13cr (see 
Fig. IU for the burst-cleaned emission confined 
in a time window of only 52 s starting at the 
onset of the first burst for the 4-50 PHA range 
(-2-20 keV) using all detector layers of PCU-2, 
the only detector operational during the observa- 
tion. This proves that the enhanced emission is 
originating from IE 1547.0-5408. Next, we revis- 
ited the enhanced emission period observed earlier 
by the PCA on MJD 54920, but in this case we did 
not find any evidence for highly increased pulsed 
emission. 




Fig. 4.— A 12 bin pulse profile of IE 1547.0-5408 
for events with PHA values in the range 4-50 (~ 2- 
20 keV) from any PCU-2 detector layer gathered 
during a 52 s time period starting at the onset 
of first burst of the Jan. 11, 2010 mini-outburst. 
Two cycles are shown for clarity. A sine fit (dashed 
line) is superposed to guide the eye and to demon- 
strate that the pulsed emission is highly symmet- 
ric. A 13cr signal is detected during this very short 
time window indicating the presence of strong pul- 
sating tail emission. 



Sim ilar behaviour has been detected bv lKaneko et al 
( 2010h using Fermi GBM data during a period of 
enhanced persistent emission of — 150 s dura- 
tion starting only 70 s after the first GBM trig- 
ger a t the onset of the Jan. 22, 2009 outburstQ 
Also, iMereghetti et al. (2009) reported on pulsat- 
ing tails lasting several seconds from two bright 
bursts occuring on Jan. 22, 2009 detected above 
80 keV by the SPI ACS aboard the INTEGRAL 
satellite. 

It is worth mentioning that given the sparse 
RXTE sampling (about 5 ks of observation time 
per week, i.e. 0.8% of the time, is devoted to ob- 
serve IE 1547.0-5408) we very likely miss a lot of 
such events i.e. periods of enhanced persistent 
emission possibly preceded by a burst event. If 
such events are accompanied by (micro) glitches, 



1 Note, that the Fermi GBM (all-sky monitor) detected the 
onset of the Jan. 22, 2009 outburst abo ut 39 minutes earlier 
than Swift BAT l|Gronwall et alj|2009l . only 4 of the sky is 



visible by the BAT at any instant) 



S 



Table 3: Incoherent timing models for IE 1547.0-5408 as derived from RXTE PCA monitoring data. 



Entry 


Start 


End 


to, Epoch 


V 


V 


V 


'v 


# 


[MJD] 


[MJD] 


[MJD,TDB] 


[Hz] 


xl0~ 12 Hz/s 


xlO" 18 Hz/s 2 


xlO" 25 Hz/s 3 


1 


54743 


54845 


54780.0 


0.482728(2) 


-22.3(11) 


-2.5(3) 


15(4) 


2 


54896 


55121 


54995.0 


0.482496(2) 


-11.0(3) 


-0.44(8) 





3 


55120 


55404 


55270.0 


0.482201(2) 


-18.1(5) 


-0.50(9) 


1.2(3) 



the (very) noisy timing behaviour (see also Sect. 
13.2.21) . inhibiting the construction of phase coher- 
ent timing models over time stretches longer than 
about 30 days, can easily be explained. 

3.2. Evolution of the pulse frequency: in- 
coherent measurements 

In each RXTE sub-observation of IE 1547.0- 
5408 a coherent pulsed sign al at a rate of ~ 0.482 
Hz (see ICamilo et al.l 120071 ) could be detected in 
the barycentered time se ries. We used a Z\- 
test (jBuccheri et al.lll983l ) search in a small, typ- 
ically 5 independent Fourier steptH wide, window 
around the predicted pulse frequency. The re- 
stricted search yielded for each sub-observation a 
best estimate of the rotation rate at the gravity 
point of the sub-observation. Because the Z 2 -test 
is distributed as a x 2 for 2 x 1 degrees of freedom 
a ler error estimate on this optimum value can 
easily be derived by determining the intersection 
points of the measured Z 2 -test distribution near 
the optimum with the value of ^ 2 max — 2.296. 

The pulse frequency measurements for IE 1547.0- 
5408 from Oct. 3, 2008 till December 25, 2010 are 
shown as bold data points along with their la 
errors in the middle panel of Fig. [TJ The mea- 
surements are relative to a coherent model (see 
later in Sect. I3.2.2[) . which could be constructed 
for the time period MJD 54855-54884 (entry #3 
in Tabled]), just after the onset of the Jan. 22, 
2009 outburst, using both RXTE PCA and Swift 
XRT data. 



2 a Fourier step is denned by Az/jpg = 1/t, in which t rep- 
resents the time span of the data period 



3.2.1. Incoherent models 

It is clear from the behaviour of the pulse fre- 
quency displayed in Fig. Q] (middle panel) that 
we can identify three time periods in which the 
spin behaviour can be described with only a lim- 
ited number of parameters. Apparent discontinu- 
ities in the timing behaviour (glitches) are present 
near MJD 54853 i.e. the onset of the Jan. 22, 2009 
outburst, and near MJD 55121 (Oct. 17, 2009). 

The frequency behaviour in each of the three 
time intervals was determined by fitting the mea- 
sured pulse frequencies and their la uncertainties 
by a simple timing model (Taylor series) with ei- 
ther 3 or 4 timing parameters (v, v, v[, 'i 1 ']): 



z {t) = v+v-(t-to) + h)-{t- 



■to) 3 
(2) 

The three incoherent timing solutions Vi nc {t) are 
listed in Table [3] and shown in Fig. Q] (middle 
panel) as solid lines together with the la uncer- 
tainty band (dashed lines), all relative to the co- 
herent timing model #3 (see Table HJ. 

3.2.2. Coherent models: phase- connected 
ephemerides 

Very accurate timing models (ephemerides) 
taking into account every rotation of the neu- 
tron star can be obtained by employing phase- 
coherent timing techniques. This requires a set 
of pulse arrival times (ToA's) which results from 
a correlation analysis of an instantaneous profile 
with a template correlation profile. For the phase- 
coherent timing techn iques used in this work w e 
refer to Section 4.1 of Kuiper fc Hermsen ( 2009f h 
In this study the correlation template is based on 
RXTE PCA events with measured energies < 4.2 
keV, corresponding to the 4-10 PHA range, be- 



9 



Table 4: Phase- coherent ephemerides for IE 1547.0-5408 as derived from RXTE PCA and Swift XRT data. 



Entry 


Start 


End 


to, Epoch 


// 


V 


V 


$o Validity range 


# 


[MJD] 


[MJD] 


[MJD,TDB] 


[Hz] 


xlO" 12 Hz/s 


xlO" 18 Hz/s 2 




(days) 


1 


54743 


54779 


54743.0 


0.48277818(3) 


-6.48(5) 


-6.46(4) 


0.4613 


37 


2 


54819 


54845 


54819.0 


0.48264979(4) 


-18.67(4) 


0.0 (fixed) 


0.7382 


27 


3 


54855 


54884 


54856.0 


0.48259525(1) 


-5.12(2) 


0.0 (fixed) 


0.0710 


30 


4 


54855 


54890 


54856.0 


0.48259518(3) 


-4.91(5) 


-0.23(4) 


0.0711 


36 


5 


54977 


55012 


54977.0 


0.48251286(6) 


-11.05(9) 


-0.18(7) 


0.5144 


36 


G 


55213 


55243 


55229.0 


0.48226226(3) 


-15.91(2) 


-1.07(12) 


0.1590 


31 


a $o is 


the phase offset to be applied to obtain consistent phase-alignment (see Eq. |4] in 


Sect. |3.2.2p 





cause in this range the morphology of the pulse 
profile is least (almost not) affected by transient 
morphology changes (see Sect]?]), which are promi- 
nently present in other bands at higher energies > 
4.2 keV. It also facilitates the inclusion of ToA's 
based on Swift XRT observations, because then 
the PCA and XRT effective energy bands are fully 
overlapping. The PCA correlation template (PHA 
4-10; < 4.2 keV) is shown in Fig. [5] adopting 15 
bins, and it is based on RXTE observations per- 
formed a couple of days after the Jan. 22, 2009 
outburst. 

Reliable phase-coherent timing models are 
given in Table |U During a 37-days period after 
the onset of the Oct. 2008 outburst IE 1547.0- 
5408 was intensively monitored by both RXTE 
and Swift (from MJD 54742-54764, only). For this 
period entry #1 of Table 0] specifies the model pa- 
rameters based on the combined PCA/XRT ToA's 
set. Thes e are c onsist ent with tho s e pub lished by 
Ng et all (|201l[ ) and llsrael et al.l (|2010l) . Note, 
the large negative value for is, indicating a rapidly 
increasing spin-down rate. 

Entries #3 and #4 of Table 0] show the tim- 
ing parameters for a period of intensive Swift and 
RXTE monitoring (just) after the second outburst 
in Jan. 2009, and are valid for a 30-days (2 pa- 
rameters) and 36-days (3 parameters) time pe- 
riod, respectively. The timi ng models are co nsis- 
tent with those r eport ed in iNg et al.1 (1201 lh and 
Bernardini et al.l ( 201 lh ■ Comparing these models 
with the one (entry #2) valid for the period just 



1.90x10* - 




Fig. 5.— The RXTE PCA pulse profile (15 bins) 
for the PHA range 4-10 (measured energies < 
4.2 keV) used as template in the Time-of-Arrival 
(TOA) correlation analysis. 

before the Jan. 2009 outburst i.e. MJD 54819- 
5484E0, indicates that a strong glitch, mainly in i>, 
occurred at a time consistent with the onset of the 
Jan. 2009 outburst. More details on this apparent 
glitch are shown later in Sect. 13.3.11 Finally, the 



3 Two Swift XRT observations in WT mode have been per- 
formed on Jan. 4 and 13, 2009, respectively. The pulse 
profile resulting from the combination of these XRT obser- 
vations adopting the timing model (entry #2 of Table |4j 
based on PCA data only, shows the expected alignment. 
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validity of coherent timing models #5 and #6 of 
Table based on PCA data only, have been ver- 
ified using independent Swift XRT data taken in 
WT-mode. Folding these Swift data yielded con- 
sistent pulse alignment. 

Also, for other time intervals with interval 
lengths of typically 30-35 days we could generate 
phase-coherent models from PCA data, mostly us- 
ing 2-d (v, v) optimization schemes applied for a 
2a range around predicted (v, v) values estimated 
from an incoherent timing model (see Table [3]). 
However, often due to the lack of independent con- 
temporaneous data from other X-ray instruments 
e.g. Swift XRlQ and given the sparse sampling 
of the RXTE observations, we can not securely 
identify the found "best" model as the true un- 
derlying rotation model. These models, not listed 
in Table |H are still useful for combining RXTE 
observations in spectral studies covering several 
sub-time intervals. 

Finally, for this source phase coherence is typi- 
cally lost after a ~ 30 — 35 day time period due to 
the sparse sampling of both RXTE and Swift ob- 
serv ations and the high level of timing noise (see 



4x1 -5 - 



e.g. ICamilo et al1l2008l) 



3.3. Glitches 

3.3.1. The characteristics of the glitch associated 
with the Jan. 22, 2009 outburst 

In order to study the Jan. 2009 glitch in detail 
we adopted standard glitch fitting techniques as 
employed e.g. in the pulsar timing software pack- 
age tempo2 vrs l.£@. We use the following, com- 
monly used, expression for the evolution of the 
frequency crossing a glitch: 

v{t) = v Q (t) + Av p + Av d e' {t " ts)/Td + Ai>(t - t g ), 

(3) 

where v^it) is the frequency of the pulsar prior 
to the glitch occuring at time t g , Av p is the fre- 
quency jump, that is permanent, while Av d is the 
part that (exponently) decays at time scale T4, and 
finally Av specifies the jump in the frequency time 
derivative. The initial (at time t — t g ) frequency 
jump Av is given by Av v + Ava- 




4 Most of the time the Swift XRT operated in PC mode, 

inadequate for timing studies of IE 1547.0-5408. 
5 See http:/ /www. atnf.csiro.au/research/pulsar/tempo2 



54820 54830 54840 54850 54860 54870 54880 
T [MJD] 

Fig. 6. — A zoom-in of the rotation frequency 
behaviour around the Jan. 2009 outburst. Fre- 
quency measurements are given by the data 
points, while the solid lines refer to coherent tim- 
ing models #2 and #3 of Table [4] All frequency 
values are relative to predictions by model #3. An 
apparent glitch occurred somewhere between MJD 
54845.02 and 54853.91, very likely simultaneously 
with the onset of the Jan. 2009 outburst at MJD 
54853.035 (TDB). 



The frequency behaviour prior and after the 
glitch (see Fig(5J) does not show any recovery term, 
and therefore we abandon the decaying part of 
the frequency jump. Assuming, that the glitch 
occurred at the time of the onset of the Jan. 
2009 outburst i.e. MJD 54853.035 (TDB) we de- 
rived the following values for the relative frequency 
jump Av/v of (1.9 ± 1.6) • 10 -6 , and frequency 
derivative jump Av/v of —0.69 ±0.07. Note, that 
the value for the frequency jump is not significant 
at face value, in contrast to the frequency deriva- 
tive jump value of Av = +(1.30 ± 0.14) • 10" 11 
Hz/s, indicating a dramatic decrease of the spin- 
down rate. Such an instantaneous decrease in 
spin-down rate, from (— 18.8±0.6) • 10 -12 Hz/s to 
(-5.12 ± 0.02) ■ 10~ 12 Hz/s, has not been seen be- 
fore from any magnetar (see the discussion section 
for a comparison with earlier reports on magnetar 
timing, Sect. [7]). 

The arrival time residuals for the best fitting 
timing model crossing the glitch, thus including 
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the frequency and frequency derivative jumps, are 
shown for the range MJD 54819-54857 in FigEl 
and indicate a proper description of the measured 
ToA's by the glitch-crossing timing model. 



Jan. '09 
Outburst 



54810 54820 54830 54840 54850 54860 

T [MJD; TDB] 

Fig. 7. — Arrival time residuals (in ms) across 
the range MJD 54819-54857 for a fit based on a 
timing model crossing the glitch epoch, that was 
hold fixed at the onset of the Jan. 2009 outburst 
at MJD 54853.035 (TDB). 



3.3.2. October 2009 glitch 

From the middle panel of FigJU showing the fre- 
quency measurements with superposed the three 
incoherent timing models, another discontinuity in 
the timing behaviour of IE 1547.0-5408 can be dis- 
cerned near MJD 55121 (Oct. 17, 2009). Unfor- 
tunately, the RXTE (and Swift) observation sam- 
pling is very sparse, such that a detailed study is 
impossible. From the incoherent models (plus un- 
certainties; see Table [3]) near MJD 55121 we can 
determine that also for this glitch a strong A/>/z> 
jump occurred of size (—0.89 ± 0.19), similar to 
the Jan. 22, 2009 glitch. More surprisingly, in this 
case the glitch is not accompanied by a radiative 
outburst (cf. top and middle panels of Fig|T]). 

3.4. Pulsed flux evolution at X-rays below 
- 10 keV 

Using either the verified (see Table [4]) or non- 
verified (see Sect l3.2.2]) coherent timing models 
we phase-folded barycentered PCA event arrival 



times according to 

*(t) = v(t-t ) + -v-(t-t ) 2 + ^i>-(t-t ) 3 -$o, 

(4) 

to obtain pulse-phase distributions, (v, i/, v) rep- 
resent the pulse frequency, first time derivative of 
the frequency and second time derivative of the 
frequency, respectively, while to is the epoch of 
the ephemeris. Consistent X-ray phase alignment 
is obtained by subtracting $o as shown in Eq. |4j 

Next, for every PCA sub-observation we fitted 
the pulse-phase distribution, based on events from 
PCU-2 with PHA values in the range [4,27] (mea- 
sured energy < 10 keV; all Xenon detector layers 
used), with a model consisting of a constant plus 
a truncated Fourier series (the fundamental and 
2 harmonics). This yielded for every PCA sub- 
observation the number of pulsed excess counts 
which in turn was converted to a (pulsed) count 
rate using the known exposure time. These PCU- 
2 (pulsed) count rates for the [4,27] PHA range 
are shown in the bottom panel of Fig. [1] for the 
time window MJD 54743-55555 (Oct. 3, 2008 - 
December 25, 2010). 

If we compare the total and pulsed X-ray flux 
measurements (both -2-10 keV) by RXTE PCA 
(cf. top and bottom panels of FigfT]), then a 
completely different evolution is shown, indicat- 
ing that the pulsed component is (at most) weakly 
related to the unpulsed or DC component. Most 
strikingly, the radiative outburst of the pulsed flux 
in October 2008 is stronger than that in January 
2009, while the total flux outburst in January 2009 
is much more intense than that in October 2008. 

Concentrating on the pulsed-flux evolution we 
notice that, i) in October 2008 the pulsed flux 
maximum is reached with a delay of about 10 
days with respect to the onset of the Oct. 3, 2008 
outburst, which is c onsistent with the findin gs by 
Ng et all (|201lh and lScholz fc Kaspil l|201ll) . and 



ii) after the start of the Jan. 22, 2009 outburst the 
pulsed flux levels appeared to vary by a factor —2 
for about two weeks, then more or less stabilize, 
and a period of gradual decrease sets in lasting 
from MJD 54868 up to MJD 55555. Assuming a 
linear decay over the full period, we find a good re- 
duced \l of 66.11/79 = 0.84 with a best fit decay 
rate of (—6.15 ±0.66) x 10~ 4 c/(s day), estimated 
for the — 2 — 10 keV band of PCU-2, which means 
a drop of about 54% over the full 688 day period 
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Table 5: Definitions of the RXTE observation time segments. 



Segment 


Begin/End Date 


Start/End 


Duration 


PCU-2 / 0-4 


nr. 


(yyyy/mm/dd) 


(MJD) 


(days) 


exposure (ks) 


1 


2008-10-03 / 2008-10-29 


54743 - 


54768 


26 


64.007 / 116.874 


2 


2008-11-01 / 2009-01-14 


54771 - 


54845 


75 


68.783 / 120.092 


3 


2009-01-22 / 2009-02-01 


54853 - 


54863 


11 


84.094 / 130.199 


4 


2009-02-06 / 2009-05-03 


54868 - 


54954 


87 


79.662 / 149.528 


5 


2009-05-11 / 2009-06-22 


54962 - 


55004 


43 


40.614 / 73.717 


6 


2009-07-14 / 2009-10-08 


55026 - 


55112 


87 


45.331 / 77.139 


7 


2009-10-16 / 2009-12-26 


55120 - 


55191 


72 


47.714 / 81.113 


8 


2010-01-17 / 2010-03-09 


55213 - 


55264 


52 


34.472 / 62.683 


9 


2010-03-17 / 2010-05-12 


55272 - 


55328 


57 


38.368 / 71.177 


10 


2010-05-17 / 2010-07-27 


55333 - 


55404 


72 


45.784 / 83.488 


11 


2010-08-04 / 2010-12-25 


55412 - 


55555 


144 


109.704 / 193.350 



PCU-2 and PCU 0-4 exposures are screened adopting default selection criteria, 
removing short duration bursts and taking into account detector break-downs. 



(see dashed line in the bottom panel of Fig[T) . 

4. Pulse profile morphology evolution at 
X-rays below ~ 30 keV 

In order to study in detail the evolution of the 
pulse profile morphology in various X-ray bands, 
we combined the data from several RXTE/PCA 
sub-observations into 11 time segments spread 
over the time period from Oct. 3, 2008 up to Dec. 
25, 2010 (MJD 54743-55555). Details of these 
RXTE observation time segments can be found 
in Table [5] The first two time segments refer to 
RXTE observations performed between the Oct. 
3, 2008 and Jan. 22, 2009 outbursts, while the 
other segments have start times beyond the Jan. 
22, 2009 outburst (see middle panel of Fig.p. 

After folding each barycentered PCA event 
time stamp according to Eq. |4] from burst and 
detector break-down free periods we generated for 
each of the time segments event matrices with 
(60 x 256) elements by binning the pulse-phase 
range [0,1] into 60 phase bins for all 256 PHA 
channels. In this process we used data from all 
three xenon layers of each PCU, which consid- 
erably improves the signal-to-noise ratio for en- 
ergies above ~ 10 keV. This allows us to better 
characterize the hard X-ray (> 10 keV) proper- 
ties. In Fig.lHthe pulse profiles of IE 1547.0-5408 
are shown for all 11 time segments adopting three 
different X-ray bands: < 4.14 keV (soft band, 
PHA 4-10), 4.14-11.22 keV (medium band, PHA 
11-27) and 11.22 - 33.82 keV (hard band, PHA 



28-80). Drastic morphology changes as a function 
of time are shown for the hard and medium energy 
bands, especially for the transitions from segment 
2 to 3, and 3 to 4. Apparently, the dramatic 
changes of physical conditions after the Jan. 22, 
2009 event (the timing glitch and radiative out- 
burst occur just prior to the start of segment-3) 
are responsible for these phenomena 

Noteworthy is that the hard X-ray bands of 
segment 1 and 2 already show pulsed emissions 
at the 4.7cr and 2.9a significance levels (apply- 
ing -tests), respectively. Therefore, the period 
just after the Oct. 3, 2008 outburst represents the 
first time period in which from this source pulsed 
emission has been detected significantly at ener- 
gies above ~ 10 keV. 

From Fig. [8] it is also clear that the soft X-ray 
pulse profiles (< 4.14 keV) exhibit least variability 
and remain relatively stable. This stability was 
the main reason for choosing the soft X-ray band 
for the To A correlation analysis (see Sect. 13 . 2 .2[) . 

4.1. Evolution of a transient hard X-ray 
pulse above ~ 10 keV 

In the hard X-ray band of segment-4 (see Fig. 
[8] right panel labeled 4) a completely new, rela- 
tively sharp, emission feature pops up near phase 
0.6 after a short transition phase during segment- 
3, which lasts only 11 days. This feature gradu- 
ally fades and is not detectable anymore beyond 
segment-7. To obtain the best statistics for this 
new component in the hard X-ray band we stacked 
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Fig. 9. — Count-rate evolution since the Jan. 22, 
2009 outburst (shown as vertical grey line) of the 
transient hard X-ray pulse (triangles) for energies 
- 11.2 - 33.8 keV (see Sect. BTTjl , The selected 
PCA time segments are indicated near the top of 
the figure. The hard pulsed component reaches 
its maximum in segment 4, which starts 15 days 
after the onset of the Jan. 22, 2009 outburst and 
lasts 87 days. Next, this component fades to unde- 
tectable levels within ~ 350 days from the onset. 
The fading trend of this component is also clearly 
visible at lower energies (~ 4.1 — 11.2 keV; dots). 
All errors are statistical only (la). 



all data from segment 4 up to and including 7, 
and fitted the resulting pulse phase distribution 
(in 60 bins) with a model composed of a Gaus- 
sian (3 free parameters) and a (flat) background 
(1 free parameter). The reduced xl of the best fit 
is 68.21/(60 - 4) = 1.22 indicating that the model 
provides an adequate description of the data. The 
resulting centroid and width (FWHM=2.354cr) of 
the best fit Gaussian are 0.587(4) and 0.23(1), re- 
spectively. 

In order to follow the evolution of the fad- 
ing hard X-ray pulse we fitted the hard X-ray 
pulse profiles of segments 3 to 11 in terms of this 
"best fit" Gaussian (fixed position and width) with 
free scale and a constant. From the fits we de- 
rived the number of excess counts associated with 



the Gaussian component, which together with the 
(screened) exposure time for the combination of 
PCU 0-4 (see 8 th column of Table EJ can be trans- 
lated to pulsed count rates. 

The quality of the nine fits was good/acceptable 
except for segments 3 and 4, yielding for both re- 
duced xl values of ~ 1.8 for 60 — 2 = 58 degrees 
of freedom. For segment-4, a slight shift to the 
left of the Gaussian template provides an excel- 
lent fit and a consistent count rate. In the case 
of segment-3 the observed profile is deviating too 
much from a Gaussian shape. A fit with a trun- 
cated Fourier series (fundamental plus 2 harmon- 
ics) to determine the unpulsed level yielded a good 
fit (reduced xl = 59.54/(60-7) = 1.12). Now, the 
pulsed excess counts derived in a phase window 
[0.426-0.748], centered on 0.587 containing a 90% 
fraction of the "best fit" Gaussian, above this dif- 
ferently derived unpulsed level resulted in a pulsed 
count rate ~ 1.6 times higher than that estimated 
by the Gaussian extraction method. Segment-3 
seems to represent a transition period starting af- 
ter the glitch epoch till a configuration/geometry 
is reached with the relatively narrow hard X-ray 
pulse at phase ~ 0.59. In Fig. [HI the hard X-ray 
pulsed count rate, adopting the Gaussian extrac- 
tion method, is shown versus time since the Jan. 
22, 2009 outburst. It is clear, that the hard X- 
ray pulsed component fades to undetectable levels 
within about 350 days from the onset of the Jan. 
22, 2009 outburst. 

Also, in spite of the difficulties in reliably es- 
timating the hard X-ray pulsed count rate for 
segment-3, we can securely state that the maxi- 
mum emission of the hard X-ray pulsed component 
is delayed with respect to the Jan. 22, 2009 out- 
burst. This is supported by non- uniformity signif- 
icance estimations of the pulse phase distributions 
in the hard X-ray band using Zf statistics for seg- 
ment 3 and 4. Although both segments have com- 
parable exposure times (see Table O , we detect 
for segment-3 only a 9.2<r pulsed signal, while for 
segment-4 we find a 27.5<r signal (cf. Fig. [8] right 
panels 3,4). 

4.2. Modelling the pulse profiles below 10 
keV 

The evolution of the pulse shape in the middle 
panel of Fig.[8](4.14-11.22 keV band) suggests that 
the profiles of segments 4 and beyond are com- 
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posed of two components: 1) a relatively stable 
one with a shape compatible with that at lower 
energies (see left panel of Fig-EJ) and 2) a transient 
Gaussian-like feature centered near phase 0.6, the 
phase of the transient pulse above 10 keV. 

In order to model the evolution of the profile 
shape in the medium band from segment 4 and be- 
yond we generated now also a template for the soft 
X-ray band (< 4.14 keV) using data from the time 
periods covering segments 8-11, a year after the 
Jan. 2009 timing glitch. This soft X-ray template 
has been obtained applying a truncated Fourier se- 
ries fit, using the fundamental and two harmonics, 
of the combined < 4.14 keV pulse phase distribu- 
tions of segment 8/9/10 and 11 (see lower two left 
panels of Fig. [8]). The resulting template is super- 
posed on the low-energy pulse phase distribution 
of segment 11 (see bottom left corner of Fig. [8]). 

Next, we fitted all < 4.14 keV pulse phase dis- 
tributions in terms of a constant and the soft X-ray 
template with free normalization, and obtained 
good/reasonable (even for segment-3) fits for all 
segments except 1 and 2 (most of the fits are su- 
perposed on the distributions shown in the left 
panel of Fig.|S|). Therefore, we conclude that the 
pre- Jan. 22, 2009 outburst soft X-ray profile dif- 
fers significantly from that of the post- Jan. 22, 
2009. The post-outburst model clearly underesti- 
mates the pulsed emission near the maximum, and 
overestimates the contribution near phase 0.6 (see 
left top panel of Fig. [5]). 

Finally, the medium-energy profiles from seg- 
ments 4 and beyond were fitted with a model com- 
posed of a constant, the soft X-ray and hard X-ray 
(Gaussian) templates both with free scale. We 
obtained acceptable fits, i.e. \l ^ 1-36 for 57 
degrees of freedom, for all time segments except 
for segment-6 (xl — 1.82). From the scale fac- 
tors of the best fits and exposure times the count 
rates for both the soft X-ray and hard X-ray com- 
ponents can be determined. For the hard X-ray 
(Gaussian) component these rates are added in 
Fig.|9l and show a similar fading trend as visi- 
ble at higher X-ray energies, meaning that the 
hard Gaussian like transient pulse extends towards 
lower X-ray energies. For segment-4 38 ±2% of the 
total pulsed emission in the medium-energy band 
can be attributed to the hard X-ray Gaussian com- 
ponent, while for segment-5 and 6 these numbers 
are 24 ± 2% and 15 ± 2%, respectively. 



5. Broad-band spectral evolution: total 
and pulsed X-ray emissions up to ~ 300 
keV 

For our spectral analyses of the total and pulsed 
emissions, we addressed the broad energy range 
from ~ 1 keV up to ~ 300 keV, exploiting the 
different capabilities of instruments aboard three 
high-energy observatories: RXTE, Swift and IN- 
TEGRAL. Data from the non-imaging RXTE in- 
struments PCA and HEXTE have been used to 
obtain spectral information on the pulsed emis- 
sion over a broad energy range, ~ 2.5 — 250 keV, 
while data from the Swift XRT instrument in PC- 
mode have been analysed to derive the total emis- 
sion spectrum at soft X-rays (0.2-10 keV). Using 
INTEGRAL ISGRI data spectral information has 
been obtained for both the total and pulsed emis- 
sion of IE 1547.0-5408 over the 20-300 keV energy 
band. Finally, during the early phase of the decay 
of the Jan. 2009 outburst JEM-X1 (3-30 keV) de- 
tected IE 1547.0-5408 up - 15 keV, and through 
imaging analysis spectral information of the to- 
tal emission could even be obtained for energies 
between 10 and 15 keV, a range not/hardly acces- 
sible for other high-energy instruments with imag- 
ing capabilities. 

All our spectral results have been corrected for 
interstellar absorption adopting a Hydrogen col- 
umn density N H of 3.12 x 10 22 cm" 2 (s e e e.g . 



Gelfand fc Gaenslerl 120071 lHalpern et all [2008, 



slightly above the total Galactic Hi column den- 
sity estimated to be in the range (1.81—2.25) x 10 22 
cm -2 ), and thus correspond to unabsorbed source 
spectra. Note, that these corrections are only ef- 
fective for photon energies below ~ 5 keV, and 
thus it is not necessary to apply these to HEXTE 
and ISGRI data. 

5.1. RXTE PCA and HEXTE pulsed X- 
ray spectra 

Analogous to the procedure followed for the 
PCA to obtain event matrices (see Sect|4j we gen- 
erated for the HEXTE data for all 11 time seg- 
ments two dimensional (burst-free) event distri- 
butions sorting on pulse-phase and energy. First, 
we verified that for the combined RXTE time seg- 
ments 4-7 we also detect with HEXTE at hard X- 
rays the same Gaussian-shaped high-energy pulse 
as measured with the PCA below ~ 30 keV. Fig. 
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Fig. 10.— HEXTE pulse- profiles for the 32.1-60.1 
(10.6cr; panel a) and 60.1-156.6 (8.6er; panel b) keV 
bands combining data from RXTE time segments 
4 — 7. The Gaussian shaped hard X-ray profile ex- 
tends up to ~ 150 keV. 



[TU] shows the HEXTE pulse profiles for the energy 
bands 32.1-60.1 keV and 60.1-156.6 keV in which 
this Gaussian profile is found at the same pulse 
phase, with non-uniformity significances 10.6 and 
8.9cr, respectively. 

Next, we derived pulsed excess counts in user 
selected energy (PHA) bands for both PCA and 
HEXTE data by fitting a truncated Fourier series 
(using the fundamental plus 2 harmonics) to the 
pulse-phase distributions measured in these bands. 
The fit function minimum plus associated error 
are used to determine the number of pulsed excess 
counts i.e. the number of counts above this mini- 
mum (=unpulsed, i.e. DC plus background level) 
for each user-selected energy band. 

For the PCA we constructed time-averaged en- 
ergy response matrices for each PCU separately 
taking into account the different (screened) expo- 
sure times of the involved PCU's during the time 
period of interest. For this purpose we used the 
ftools version 6.4 programs pcarsp and addrmf. To 
convert PHA channels to measured energy val- 
ues, Epha, for PCU combined/stacked products 
we also generated a weighted PCU-combined en- 
ergy response matrix. 

For HEXTE we employed cluster A and B 
energy-response matrices separately, taking into 
account the different screened on-source exposure 
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Fig. 11. — Spectral evolution from Oct. 2008 to 
Dec. 2010 of the pulsed (unabsorbed) emission of 
IE 1547.0-5408 in the 2.5-150 keV band as mea- 
sured by RXTE PCA (filled circles) and HEXTE 
(filled squares) for time segments 1-11 (see Table 
[5] and middle panel Fig.[lJ . Power- law plus cutoff 
power-law model fits are superposed for segments 
3, 4 and 5/6/7 (see Table [6]), during which sig- 
nificant pulsed emission was detected by HEXTE 
above ^15 keV, otherwise only power-law- model 
fits are shown below 10 keV. The pulsed hard X- 
ray emission (> 10 keV) becomes maximal during 
segment-4, showing therefore a delayed increase 
with respect to the onset of the Jan. 2009 out- 
burst, before a gradual decline sets in. 



times and the reduction in efficiency in case of off- 
axis observations. The on-source exposure times 
for both clusters have been corrected for the con- 
siderable dead-time effects. 

To translate the derived pulsed excess counts 
to photon flux values we employed these response 
matrices in (forward folding) spectral fitting pro- 
cedures adopting simple underlying photon emis- 
sion models like a power-law, a combination of two 
power-laws or a power-law plus an (exponential) 
cutoff power-law. 

The results of the spectral fits to the combined 
PCA/HEXTE data setU for the pulsed emission 



'We applied the following energy-independent normaliza- 
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Table 6: Characteristics of the spectral fits, adopting a combination of a power-law and a cutoff power-law, 
to RXTE PCA and HEXTE (2.5-150 keV) pulsed flux data. 



RXTE E ki ri k 2 r 2 E c fI~ w 

segm. (kcV) xlO 6 X 10 s (keV) 

3 12.9295 1.10 ±0.11 -4.63l°;^ 1.49 ±0.32 

4 12.8760 2.68 ±0.21 — 3.92±g;g| 1.86 ±0.33 
5/6/7 12.8875 2.50 ±0.15 — 3.92±3".ol 0.43 ±0.15 

1) The (photon flux) fit model is: F 1 (E 1 ) = fci ■ (_E 7 /E ) ri + k 2 ■ (E 1 /E ) F ^ ■ exp(-£ 7 /,E c ) 

2) The normalisations k\ and k 2 are in units ph/cm 2 s keV at pivot energy Eq. 

3) The unabsorbed energy fluxes F^T 10 and F J, 0-150 are in units 10 -11 erg/cm 2 s 

4) Quoted errors are for a 68.27% confidence level (la) 



±0.37I° 2° 44.0 ± 10.4 



±0.591q 25 33.8 ±5.9 1.28 
+ L55 -o'.42 25.3 ±4.i 
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of all 11 segments along with the (reconstructed) 
flux measurements are shown in a vF v representa- 
tion in Fig. [TT] for the range ~ 2.5 — 150 keV. 
In this figure the HEXTE results for segments 

I, 2 and 11 are left out for clearity, because the 
non-detections in the different energy bands corre- 
spond to 2a upper-limit levels comparable to those 
shown for segment-8/9/10. However, HEXTE did 
detect pulsed emission (3.6c) in its integral energy 
band, ~ 15 — 250 keV, for segment-1, confirming 
the PCA findings (see Sect. gj. 

Below ~ 10 keV the pulsed spectra are all, irre- 
spective of the time segment, very soft with pho- 
ton indices T in the range —[3.9-4.2] (except for 
segment-3, T ~ —4.6). Also, a significant decrease 
by a factor of ~ 2 in normalization is shown for 
the time period between the two outbursts, mov- 
ing from segment- 1 to 2. Then, promptly after 
the Jan. 2009 outburst a strong increase in nor- 
malization by a factor of ~ 2 is observed, mov- 
ing from segment- 2 to 3, followed by a gradual 
decrease down to the level shown for segment- 

II. The whole picture for energies below ~ 10 
keV is consistent with the behaviour in the coarse 
(^2 — 10 keV) energy band shown in the bottom 
panel of FigfT] 

Above ~ 10 keV the variations are more drastic. 
Since the initial activation of pulsed high-energy 
emission after the Oct. 2008 outburst this compo- 



tion factors of 0.912 (PCA) and 1.087 (HEXTE) to the 
spectral measurements in order to obtain co nsistent high- 
energ y Crab pulsed fluxes (see Sect. 3.4 of iKuiper et al.l 



nent reached undetectable levels during segmcnt- 
2, but increased dramatically after the Jan. 2009 
outburst, reaching eventually a maximum during 
segment-4, after which a gradual decrease sets in 
resulting to undetectable levels beyond segment-7. 

During the periods with significant pulsed hard 
X-ray emission the best spectral description over 
the full 2.5-150 keV band consists of a combina- 
tion of a power-law plus a cutoff power-law model 
(see Table |6] for quantitative information) . The 
improvement relative to a model composed of two 
power-laws is about 3.5cx, considering the com- 
bination of the three periods, segments 3, 4 and 
5/6/7, during which significant pulsed hard X-ray 
emission has been detected. 

For the cutoff power-law model component, 
F(E^) = k ■ (E 1 /E ) T ■ exp(-E 1 /E c ) with k the 
normalization in ph/cm 2 s keV, Eg the pivot en- 
ergy, r the power-law (photon) index and E c the 
cutoff energy, the photon index at a specific energy 
r* := is given by V* = T - E y /E c . There- 

fore, for this model component the maximum in a 
vF v spectral representation (maximum power per 
energy decade) is found at E™ ax = (T + 2) ■ E c . 
For the spectral fits of segments 3, 4 and 5/6/7 the 
best fit parameters for T and E c yield vF v maxima 
for energies in the range 90-105 keV with an un- 
certainty of about 20 keV. We stress that our data 
donot constrain the extrapolation of the spectra, 
the spectral shape, above these energies. 

IE 1547.0-5408 exhibits in time segment 4, ex- 
tending over 87 days, the maximal pulsed flux at 
energies above 10 keV. An important parameter 
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for theoretical modelling is the time between the 
onset of the outburst/timing glitch and the epoch 
at which the maximal flux of the non-thermal tran- 
sient component is reached. We used the high 
count rate in the RXTE PCA to divide segment 4 
in three time intervals to study the spectral evolu- 
tion in more detail. The normalisation of the spec- 
tra appeared maximal in two intervals between 37 
and 101 days after the onset of the Jan. 2009 out- 
burst. We refer to this delay as 70 ±30 days. 

5.2. Swift XRT soft X-ray spectra of the 
total emission 

To study the evolution of the total (=pulsed 
plus unpulsed) emission spectrum of IE 1547.0- 
5408 at soft X-rays (0.5-10 keV) since the Jan. 
2009 outburst we used Swift XRT data obtained 
in PC-mode, allowing full 2D-imaging informa- 
tion (see Sect. l2.3p . This information is necessary 
to get rid of the emissio n from the dust scatter 
rings ( Tiengo et al.ll2010h which were prominently 
present during the early part of the decay phase 
after the Jan. 2009 outburst. 

The Swift observations selected for this analysis 
are chosen such, that these cover (sometimes ap- 
proximately) time periods for which hard X-ray 
INTEGRAL observations have been performed 
(see Tableland Sect.ESJfor the INTEGRAL and 
Swift observations, respectively), supplemented 
with some recent observations (e.g. those with 
obs. identifiers 00090404019, Sept. 28, 2010, and 
00090404027, Feb. 25, 2011) to follow the late 
time evolution. 

We selected events from a 60" circular aper- 
ture centered on IE 1547.0-5408 using (default) 
cleaned event lists. For this selection we exe- 
cuted procedure xrtmkarf (vrs. 0.5.8 embedded 
in the HEAS0FT instrument specific software anal- 
ysis environment) to obtain appropriate sensitive 
area information, stored in so-called ancillary re- 
sponse files (arf), which takes into account the 
reduction in sensitive area for off-axis observa- 
tions and the finite angular size of the source ex- 
traction region (about 90% of the source counts 
falls within a 60" circle centered on the source). 
The energy redistribution matrix (rmf) used was 
swxpc0tol2s6_20070901v011.rmf . 

We also selected a background region composed 
of an annular region with a 180" inner- and 240" 



outer radius centered on IE 1547.0-5408, 7 times 
larger than the source region. From this we es- 
timated the number of background counts in our 
source region, which was < 1% of the total number 
of counts confined in our source region, and there- 
fore no background correction has been made in 
subsequent spectral fitting procedures. 

Sorting the selected (source region) events on 
energy we fitted the resulting event distribution 
in measured energy space to a certain model pho- 
ton spectrum attenuated by interstellar absorption 
(see Sect. [5]) in a forward folding procedure us- 
ing the appropriate response information (arf and 
rmf response files). As model (source) spectrum a 
combination of a black-body (BB) and power-law 
(PI) was used which can adequately reconstruct 
(empirically) the measured spectral distributions 
below 10 keV. These model fits were also succes- 
fully applied and discuss e d for I E 1547.0- 5 408 b y 
e.g. iGelfand fc Gaenslerl d2007lh iNg et all (l201lh : 
iBernardini et all (|201lh : IScholz fc Kaspi I (|201lh . 

The spectral fit parameters adopting this com- 
posite BB+P1 model are given in Table [7] and the 
flux measurements for the 0.5-10 keV range re- 
constructed from spectral fits are shown in Fig. [TBI 
Energy fluxes provided in Table [7] are unabsorbed 
and specified for the 1-10 keV band to facilitate 
compa risons with the results of IScholz fc Kaspi I 
(|201ll ). All fits have acceptable reduced x 2 -values 



near 1 assuming a fixed Nh of 3.12 x 10 cm . 
The fit values for the BB temperature kT and 
power-law index T as a function of time since the 
start of the outburst are consistent with those 
reported by e.g . 



Scholz fc Kaspi 



.Bernardini et al 
I 201lh 



(|201lt ) and 



Table [7] also confirms a 
clear softening of the power-law index with time 
and a decreasing flux below 10 keV (also visible 
in the spectral shapes below 10 keV in Fig. [T3| . 
This hardness-intensity c orrelation has been dis - 
cussed previously, e.g. bv lScholz fc Kaspil(|201lh . 
however, our spectral fits to the combined Swift- 
INTEGRAL spectra over the broad energy band 
1-300 keV provide a different, alternative decom- 
position and hardness-intensity correlation (Sect. 
El. 



5.3. INTEGRAL soft 7-ray total emission 
spectra 

The ISGRI observations (science windows) se- 
lected for our study all have observation dates be- 
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Table 7: Characteristics of the spectral fits, adopting a combination of a power-law and a black-body, to solely 
Swift XRT-PC (0.5-10 kcV) total&xx data. Note, however, that the T values change, when the INTEGRAL 
flux values above 20 keV are included in the fits (see Table [8]) . 



Swift XRT 


Date 


MJD 


range 


kT 


r 


pl-io 

E 


pl-io 

r E 


color 


obs. 








(keV) 




(Pl-model) 


(BB-model) 


in Fig|T3l 


00030956034 


29/01/2009 


54860.005- 


-54860.343 


0.606 ± 0.007 


-1.14 ±0.04 


5.07 ±0.50 


4.48 ±0.32 


orange 


00030956039 


4/02/2009 


54866.702- 


-54866.974 


0.564 ± 0.006 


-1.32 ±0.04 


4.41 ±0.37 


4.26 ±0.31 


purple 


00030956051 


19/08/2009 


55062.069- 


-55062.150 


0.775 ±0.018 


-1.99 ±0.16 


1.71 ±0.45 


2.52 ±0.36 


dark blue 


00090404003 


12/04/2010 


55298.509- 


-55298.938 


0.728 ±0.014 


-2.84 ±0.20 


0.93 ±0.22 


1.36 ±0.16 


aqua 



1) The (photon flux) fit model is: F(E~,) = a bb ■ E^/(exp(Ej/kT) - 1) + a pl ■ 

2) The energy fluxes F^T 10 are unabsorbed (Njj = 3.12 ■ 10 22 cm -2 ) and in units 10 — 11 erg/cm 2 s for the 1-10 keV band 

3) Quoted errors are for a 68.27% confidence level (ltr) 



tween Oct. 3, 2008 and Apr. 3, 2010 and instru- 
ment pointings within 14?5 from the IE 1547.0- 
5408 sky location. This ensures that (a part 
of) the detector plane is illuminated by the tar- 
get. The resulting list is further screened on er- 
ratic count rate variations, indicative for particle- 
induced effects due to Earth-radiation-belt pas- 
sages or solar-flare activities, by inspecting visu- 
ally the count rate in 20-30 keV band versus time. 
Science windows showing erratic count-rate varia- 
tions are excluded for further analysis. Moreover, 
we ignored time intervals during which bursts 
(from any source in the field of view) occurred. 
This selection is particularly important for INTE- 
GRAL revolutions 767-772 during which copious 
numbers of bursts from IE 1547.0-5408 have been 
detected. The set of INTEGRAL revolutions used 
along with the effective on-axis exposure times is 
given in Table [2] 

The total emission spectrum at soft gamma 
rays / hard X-rays above 20 keV can be derived 
by exploiting the arcminute imaging capability of 
ISGRI. We generated skymaps for each selected 
science window (see Sect. l2.2|) in different energy 
bands covering the 20-300 keV energy window. 
From these maps sky mosaics were assembled us- 
ing dedicated imaging software tools embedded in 
the OSA 7.0 analysis environment. We chose in 
most cases a logarithmically binned energy grid 
with 10 bins covering the 20-300 keV range, but 
also 4 logarithmic energy bins over the 20-150 
keV range when we expected low flux levels for 
IE 1547.0-5408 (e.g. for INTEGRAL Revs. 899- 



912). The final product of the imaging analy- 
sis provides the count rate, its variance, exposure 
and significance for all energy bands over the (de- 
convolved) mosaiced sky field. As an example 
Fig. [12] shows the significance map (20-150 keV) 
in a wide field centered on IE 1547.0-5408 com- 
bining observations from INTEGRAL revolutions 
899-912, performed more than a year after the Jan. 
2009 outburst. IE 1547.0-5408 is clearly detected 
in this mosaic at a level of ~ 9.5<r. 

The total emission spectrum of IE 1547.0-5408 
can be derived from the count-rate and vari- 
ance maps by extracting the (dead-time cor- 
rected) rates and uncertainties at the location 
of IE 1547.0-5408. These values are normalized to 
the count rates measured for the total (nebula and 
pulsar) emission from the Crab in similar energy 
bands, using Crab calibration observations (in 
5x5 dither mode) performed during INTEGRAL 
revolution 102. From the ratios and the photon 
spectrum of the total emission from the Crab, we 
can derive the total high-energy photon spectrum 
of IE 1547.0-5408 (pulsed and any unpulsed point 
source component). For the total Crab photon 
spectrum w e used the broken-power-law spectrum 
derived by Ijourdain fc Rogues! ( 20081 ) based on 
INTEGRAL-SPI observations of the Crab at en- 
ergies between 23 and 1000 keV. The latest Crab 
cross calibrations between SPI and IBIS-ISGRI 
provided consistent results. 

Fig.[T3lshows the ISGRI flux measurements (20- 
300 keV) for different sets of observations per- 
formed between Jan. 24, 2009 and April 3, 2010 
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Table 8: Spectral properties of the total emission of IE 1547.0-5408 across the 1-300 keV range for various 
INTEGRAL observation periods (orbital revolution numbers, see Table[2]) after the timing glitch in Jan. 2009, 
derived by fitting a black-body plus power-law model to (nearly) contemporaneous INTEGRAL ISGRI and 
Swift XRT data. INTEGRAL revolutions 46-411 represent a total of 4 Ms of observations of IE 1547.0-5408 
before March 2006. 



INTEGRAL 
Obs. period 


Flux 
20-150 keV 


r 


kT 
(keV) 


BB-FluxW 
2-10 keV 


Pl-FluxW 
2-10 keV 


767 
768-772 
782-79lt 3 ) 
840-850 
899-912 
46-411' 2 ' 3 ! 


(2.52±0.37)E-10 
(1.85±0.18)E-10 
(2.15±0.14)E-10 
(0.84±0.27)E-10 
(0.80 ± 0.22)E-10 
< 0.15E-10 


-1.41 ±0.06 
-1.45 ±0.04 
-1.27±0.11 
-1.22 ±0.10 
-0.87 ±0.07 
-1.35 (fixed) 


0.621 ± 0.008 
0.566 ±0.007 

0.739 ± 0.014 
0.654 ±0.010 


(2.51 ±0.22)E-11 
(2.35±0.19)E-11 

(2.23±0.27)E-11 
(1.18±0.11)E-11 


(4.49 ± 0.23)E-11 
(3.66±0.17)E-11 

(0.91 ±0.17)E-11 
(0.35 ± 0.06)E-11 



fluxes are unabsorbed adopting Njj = 3.12 X 10 22 cm 2 and in units erg/cm 2 s 
' 2 1 Upper limit at 2a confidence 

PlOnly ISGRI data used adopting power-law model 



(Revs. 767-912; see Table [2]). For the (only) 
100 ks ToO observation performed during 8-10 
Oct., 2008 (Rev-731; a couple of days after the 
Oct. 3, 2008 outburst) no significant flux values 
can be derived with 2a upper-limits (not shown 
in Fig. [T5|) consistent with the pulsed flux detec- 
tions by RXTE PCA and HEXTE. We added in 
Fig. IT3l also the 2a upper-limits on the IE 1547.0- 
5408 flux derived from a deep 4 Ms mosaic (see 
e.g. iKuiper et al.ll2008l ) using INTEGRAL obser- 
vations performed between revolution 46 and 411 
(March 2, 2003 - Feb. 24, 2006) targeting on 
PSR J1617-5055, which is located within 5?5 from 
IE 1547.0-5408. These upper-limits are at least 10 
times lower than the flux levels reached during IN- 
TEGRAL revolutions 767-791 in which the hard 
spectral tail is detected most significantly. 

During the intensive monitoring campaign 
(Revs. 767-772) of IE 1547.0-5408 the source was 
also for - 264.8 ks within 3?75 of the JEM-X1 
pointing axis. In the combination of these obser- 
vations JEM-X1 clearly detected IE 1547.0-5408 
at a 8.6cr level between 3.04 and 15.04 keV. Flux 
values are obtained using similar procedures as 
applied for ISGRI i.e. count rates are expressed in 
Crab (total) flux units and subsequently converted 
to flux units, however, now we have to take into 
account the different column densities along the 
line of sight for the Crab (N H of 3.2 x 10 21 cm" 2 ) 
and IE 1547.0-5408 (N H of 3.12 x 10 22 cm" 2 ). 
The four JEM-X1 flux measurements (3 detec- 
tions and 1 upper-limit) are shown in Fig.[T3] as 



filled squares. 

5.4. INTEGRAL and Swift XRT total- 
emission spectral evolution 

To extend the spectral coverage down to ~ 1 
keV we combined the ISGRI(/JEM-X1) observa- 
tions with the (nearly) contemporaneous Swift 
XRT observations performed in PC-mode (see 
Sect. 15. 2[) . This allowed us to study the spectral 
shape and spectral evolution over the broad 1— 
300 keV interval. We found that a fit with a 
black-body plus (single) power-law model to the 
flux measurements over the total energy range 
provided a statistically good fit. No fudge fac- 
tors were used in these fits to absorb calibration 
uncertainties between different instruments. The 
model fits are shown in Fig.[T3] The spectra of 
the two late-time Swift-XRT observations of Sept. 
28, 2010 and Feb. 25, 2011 are not shown for 
clarity, because these are fully overlapping with 
the spectrum (aqua-coloured) of the observation 
performed on April 12, 2010 slightly after INTE- 
GRAL revolutions 899-912. 

Characteristics of these XRT-ISGRI combined 
spectral fits are listed in Tabled] The first thing 
to notice is that also for this broad energy band 
two spectral components, a BB plus a power-law 
model, a ppear to be sufficient t o fit the flux values 
(see also lBernardini et al.ll201lh . The fit values for 
the BB temperature kT (Table|8|) are close to those 
derived for solely 0.5-10 keV band Swift-XRT data 
(Table [7]), however, the power- law photon indices 
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E [MeV] 

Fig. 13. — Spectral evolution after the timing 
glitch/outburst in Jan. 2009 up to April 2010 
of the total (unabsorbed) emission of IE 1547.0- 
5408 in the 1-300 keV band as measured by 
Swift-XRT (PC-mode; filled circles), INTEGRAL 
JEM-X (filled squares) and ISGRI (filled dia- 
monds) . Black-body plus power-law-model fits for 
the combined fluxes are shown. The different col- 
ors indicate the INTEGRAL observations (revo- 
lutions) with (nearly) contemporaneous Swift ob- 
servations. Note, that for the period covering 
Revs. 782-791 no XRT observations have been 
performed in PC-mode. Black 2a upper-limits are 
shown for the combined exposure of ~ 4 Ms of all 
INTEGRAL observations (prior to the Oct. 2008 
outburst) from March 2003 to Feb. 2006. 



r in the third column of Table|S] are very differ- 
ent, ranging from ~ —1.4 after the Jan. 2009 
glitch/outburst to ~ —0.9 almost 14 months later, 
while the 20-150 keV flux is fad ing by a factor of 



~ 3. In agreement with this, Iden Hartog et al 
(|2009h . using - 200 ks of INTEGRAL (near real 
time) ToO data taken 2-8 days afte r the Ja n. 2009 
glitch/outburst, and lEnoto etal. ( 2010al) . using 
Suzaku data taken ~ 7 days after the Jan. 2009 
glitch/outburst, consistently reported a power-law 
index of 1.5 during the first half of INTE- 
GRAL observation p eriod Revs. 767-772 (see Ta- 
ble [T]). Furthermore. lEnoto et all ( 2010al Suzaku 



data) and iBernardini et al. ( 2011 . Chandra and 
INTEGRAL data) also derived for energies ~ 0.5- 
100 keV within statistical uncertainties the same 
fit result as given in Table|3]for epochs in INTE- 
GRAL Revs. 767-772. Thus, the power-law com- 
ponent resulting from the fits below 10 keV, is not 
required when the broad energy range is consid- 
ered. The "new" hard power-law component ex- 
plains the total hard X-ray emission above ~ 10 
keV and also contributes significantly to the emis- 
sion below 10 keV. Moreover, the index is hard- 
ening with decreasing flux, the opposite trend to 
what is visible i n Tab le [7] and discussed in e.g. 



Scholz fc Kaspi I (|2011l ) for spectral fits restricted 



to energies below 10 keV. 

The fading trends of the total fluxes at soft X- 
rays and hard X-rays since the Jan. 2009 outburst 
are shown in Fig. [14] The total (unabsorbed) soft 
X-ray flux is the sum of the black-body and power- 
law components (last 2 columns of Table[5]). The 
power-law contribution at soft X-rays is negligible 
after April 12, 2010. Namely, late-time Swift XRT 
measurements of the total soft X-ray spectrum of 
IE 1547.0-5408 demonstrate that a single black- 
body component is sufficient to describe the spec- 
trum adequately. The (late-time) flux measure- 
ments of Swift observations 0009040419 (Sept. 28, 
2010) and 00090404027 (Feb. 25, 2011) are added 
in Fig. [14] The corresponding (unabsorbed) 2- 
10 keV fluxes and temperatures are (1.67(9) • 
10~ n erg/cm 2 s, kT=0.671(5)) and (1.06(7) -lO -11 
erg/cm 2 s, kT=0. 600(6)), respectively. 

It is evident from Figs. [13] and [14] that the to- 
tal source fluxes are maximal during the first IN- 
TEGRAL ToO observation after the Jan. 2009 
glitch/outburst epoch (Rev. 767, starting 2 days 
after the Jan. 2009 outburst) for energies below 
10 keV as well as for the soft gamma-rays/hard 
X-rays (20-150 keV). It is remarkable that for the 
latter total non-thermal emission we do not see ev- 
idence that maximal luminosity is reached about 
70 ± 30 days after the glitch (in RXTE and IN- 
TEGRAL segments 4), a behaviour we found for 
the transient pulsed component for energies above 
~ 20 keV. Apparently, the increase of the pulsed 
emission above 20 keV is delayed with respect to 
that of the unpulsed-emission above 20 keV. The 
contribution of the pulsed component to the total 
flux (=pulsed fraction) in the 20-150 keV band 
during INTEGRAL segment 4 amounts 34 ± 9%. 
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Fig. 14. — Evolution of the total (unabsorbed) 
energy flux of IE 1547.0-5408 since the Jan. 2009 
outburst in two different energy bands: 2-10 keV 
(Swift- XRT; solid squares; right flux scale) and 20- 
150 keV (ISGRI; filled triangles; left flux scale). 
The onset of the Jan. 2009 outburst is indicated 
by a vertical line, and the INTEGRAL observation 
segments (see TableEJ are shown in the top part 
of the figure. The horizontal line represents the 
2a upper limit for the 20-150 keV flux as derived 
from INTEGRAL observations performed during 
Revs. 46-411 (before March 2006). 



At the start of the Jan. 2009 outburst this frac- 
tion was 31 ± 7%, and about 240 days later during 
INTEGRAL segment 5 45 ± 14%. 

6. Summary 

In this paper we presented new and detailed 
characteristics of the persistent (=non-burst) 
emission of IE 1547.0-5408 between ~ 1 and 300 
keV analyzing data collected over a period of 27 
months, starting with the onset of an outburst 
in October 2008, covering the epoch of the Jan- 
uary 2009 outburst, and ending two years later 
in January 2011. Detailed results are derived on 
the timing characteristics, the evolutions of the 
total and pulsed persistent emissions, and pulse 
profiles. In particular, the evolution of the total 
and pulsed emission of a magnetar for energies 
above ~ 10 keV after the event of a major timing 



glitch accompanied with a radiative outburst has 
been revealed for the first time. We used RXTE 
PCA and Swift XRT monitoring data to derive 
incoherent and coherent timing models for the 
pulse frequency evolution between October 2008 
and January 2011. RXTE PCA and HEXTE data 
were used to derive the evolution of the pulse pro- 
file, flux and spectrum of the pulsed emission from 
~ 2.5 to 150 keV. Finally, we exploited the imag- 
ing capabilities of Swift XRT and INTEGRAL 
ISGRI to derive the evolution of the flux and 
spectra of the total emission for energies from ~ 1 
to 300 keV after the January 2009 outburst. What 
follows is a summary of the main results. 

6.1. Total flux evolution of IE 1547.0-5408 
below 10 keV 

1) Using RXTE-PCA data we derived the to- 
tal flux evolution for energies — 2 — 10 keV from 
Oct. 3, 2008 up to Dec. 25, 2010 showing the 
outburst in October 2008 and the much more in- 
tense outburst in January 2009. The flux decay 
of the latter has still not reached its quiescent 
level ~2 years after the out burst (Sect.l3.ll up - 
per panel of Fig [Tj). Se e also Israel et afl ~i 2010h ; 
Bernardini et al.l (jioifh ; IScholz fc Kaspi I (|201lh . 

2) We can satisfactorily describe the decay of 
the total flux after the intense Jan. 22, 2009 
(MJD 54853.911) outburst with a rapidly decay- 
ing power-law component, dominating only dur- 
ing a period of ~5 days after the outburst, and 
a slowly decaying exponential component with a 
typical time scale of —200 days (Sect. [3~7\ . 

3) In the tail of the intense Jan. 22, 2009 out- 
burst we detected two mini-outbursts, at MJD 
54920 and at MJD 55207, with enhanced emis- 
sion over periods of —450 s and —130 s, respec- 
tively, and superposed with short duration bursts. 
In the tail of the latter outburst we found in a 
time window of only 52 s a pulsed signal (energies 
-2 - 20 keV) of -13 a significance (Sect JRTFty . 
The latter event might be similar to the persis- 
tent, pulsed e missio n enhancement reported by 
Kaneko et al. ( 2010h using Fermi GBM for a time 
interval just after the onset of the Jan. 2009 out- 
burst of IE 1547.0-5408. 
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6.2. Evolution of the pulse frequency of 
IE 1547.0-5408 and glitch detections 

1) From a timing analysis of RXTE PC A sub- 
observations we derived the evolution of the pulse 
frequency from Oct. 3, 2008 till Dec. 25, 2010, 
based on incoherent measurements (Sect. . HOI mid- 
dle panel of Fig.\l§. 

2) Two apparent discontinuities in the pulse- 
frequency evolution are visible, near MJD 54853 
(the onset of the Jan. 22, 2009 outburst) and near 
MJD 55121 (Oct. 17, 2009). For three time in- 
tervals, before, in between and after the epochs of 
these timing discontinuities, we derived simple in- 
coherent timing models (Taylor series) with either 
3 or 4 timing parameters (i/, i>, v'\) ( Sect. [3~2~7\ 
Tabled- 

3) Phase- coherent timing models have been 
constructed for 6 time intervals, each of duration 
-30 - 35 days, using RXTE PCA and Swift XRT 
overlapping observations. Most importantly, one 
phase-connected ephemeris was derived after the 
onset of th e Oct. 2008 outbur st (see also lNg et al 



to derive the pulsed flux evolution for energies 
-2-10 keV from Oct. 3, 2008 up to Dec. 
25, 2010 (Sect.W^l an d bottom panel of Fig.\$. 

W fee 



(1201 ll ) and llsrael etHI (|201dh . and one just before 
the January 2009 outburst and two directly after 
this second outburst (Sect. \3.2.B and Table^. 

4) Applying standard glitch fitting techniques 
we measured at the time of the onset of the Jan. 
2009 outburst (i.e. MJD 54853.035) a frequency 
jump Av/v of (1.9 ± 1.6) • 10 -6 , and frequency 
derivative jump Av/v of —0.69 ± 0.07. The fre- 
quency derivative jump value of Av = +(1.30 ± 
0.14) • 10~ n Hz/s indicates a dramatic sudden de- 
crease of the spin-down rate, not seen before from 
any magnetar ( Sect. [3~3~l\ see also discussion sec- 
tion, Sect. [7]). 

5) For the discontinuity in the timing behaviour 
of IE 1547.0-5408 near MJD 55121 (Oct. 17, 
2009), we estimated from the incoherent timing 
models that a strong Avjv jump occurred of size 
(-0.89 ±0.19), similar to the Jan. 22, 2009 glitch. 
Note, this glitch was not accompanied by a radia- 
tive outburst (Sect. \3.3.2\ . 

6.3. Pulsed flux evolution of IE 1547.0- 
5408 below 10 keV 

1) Applying our derived timing models for 
phase-folding barycentered RXTE PCA event ar- 
rival times, we derived pulse-phase distributions 
for every PCA sub-observation. This allowed us 



Scholz fc Kaspi I ([20111 ) showed a consistent evo- 
lution covering the first — 150 days, including the 
two outburst epochs. 

2) The evolution of the pulsed flux (energies 
— 2 — 10 keV) behaves very different compared 
to that of the total flux below 10 keV: the radia- 
tive outburst of the pulsed flux in October 2008 
is stronger than that in January 2009, contrary to 
the findings for the total flux (Sect. \3~J\ and Fig. [fj. 

3) In October 2008 the pulsed flux (-2-10 
keV) reaches its maximum — 10 days after the 
start of the outburst of the to tal emissi o n (bo t- 
tom panel of Jfo.[H). See also iNg et all ( 2011 ); 



Scholz fc Kaspi I poill) . 

4) After the start of the January 2009 outburst 
the pulsed flux (—2 — 10 keV) was unstable, vary- 
ing by a factor of — 2 for two weeks and then a 
linear decay set in over a period of almost 700 days 
with a decay rate of (— 6.15±0.66) x 10~ 4 c/ (s day) 



(Sect. \3~^\ and bottom panel of Fig.{^. 

6.4. Evolution of the pulse-profile mor- 
phology below — 30 keV 

1) Using RXTE PCA monitoring data we con- 
structed for 11 time segments between Oct. 3, 
2008 and Dec. 25, 2010 pulse profiles in three X- 
ray bands: < 4.14 keV (soft band), 4.14 - 11.22 
keV (medium band) and 11.22 - 33.82 keV (hard 
band) (Seci.g] and Fig^. 

2) Between the onset of the October 2008 out- 
burst and that of January 2009, in the soft and 
medium bands a single broad pulse is detected 
with a peak at pulse phase — 0.9. In the same time 
window first evidence was found for weak pulsed 
emission in the hard band above 11 keV (Sectj^] 
and Fig.[fJ§. 

3) Drastic changes in pulse shape are found as 
a function of energy and of time after the glitch 
epoch and start of the radiative outburst in Jan- 
uary 2009: most notably, a transient hard X-ray 
pulse appears in the phase distribution after — 
11 days, that can be well described with a Gaus- 
sian shape, centered at phase — 0.59 and width 
(FWHM) 0.23. In the hard band above - 11 keV 
this high-energy pulse is the single pulse detected; 
in the medium band this high-energy pulse is seen 
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in addition to the broad pulse detected in the soft 
band, which is similar but not identical to that 
detected in the soft and medium bands before the 
January 2009 outburst (Sect.\4~J\ FigM Sect.\J^\ 
Fig.\^. 

4) The new high-energy pulse reaches maximal 
flux in a time interval between 15 and 102 days 
after the onset of the January 2009 outburst and 
smoothly fades to undetectable levels within ~ 350 
days from the onset {Fig.\^. 

5) The time interval of ~ 11 days after the 
glitch in January 2009 appears to be a transition 
period in which the pulse profile changes from the 
single broad soft pulse before the glitch to the 
structured pulse consisting of a soft and a tran- 
sient hard component with peak phase separation 
~ 0.3. This ~ 11-day time interval is also the 
period in which the pulsed flux (~ 2 — 10 keV) 
was unstable, varying by a factor of ~ 2 before 
a linear decay set in over a period of almost 700 
days (Fig.\^ and for variability Sect. \ 3.4] and bot- 
tom panel of Fig.[I$) . 

6.5. Evolution of the pulsed spectra for X- 
ray energies 2.5 - 150 keV 

1) Exploiting the simultaneous monitoring with 
the PCA and HEXTE aboard RXTE we derived 
the spectra of the pulsed emission over the broad 
energy range 2.5 - 150 keV in 11 time segments 
between Oct. 3, 2008 and Dec. 25, 2010 (Sect.UTl] 
and Fia. [T7\i . 

2) Below ~ 10 keV all pulsed-emission spectra 
are very soft with photon indices T in the range 
-[3.9 - 4.2], except in the first time segment (11 
days; RXTE segment-3) after the Jan. 2009 out- 
burst (r ~ —4.6). In the time period between the 
two outbursts, the flux (normalization) was high- 
est during the first segment (26 days), decreased 
by a factor of ~ 2 in the second segment (75 days) 
before the outburst in Jan. 2009. Promptly af- 
ter the latter outburst, over a period of 11 days, 
the normalization was again higher by a factor ~ 
2, followed by a gradual decrease by a factor ~ 2 
over almost 700 days (Sect. \5. 1\ and Fia. [Tl\i . 

3) Around ~ 10 keV the pulsed spectra change 
from very soft spectra to hard non-thermal spec- 
tra, with jumps in spectral photon index in the 
range 2.5 - 3.0 (Fig.\ similar to what has been 
reported for the persistent pulsed emission of some 



AXPs (e.g. lKuiper et al.ll2006l) . 

4) Above ~ 10 keV the variations are more 
drastic than at lower energies. Between the Oct. 
2008 and Jan. 2009 outbursts the pulsed emission 
above 10 keV was too weak to derive a spectrum. 
Directly after the Jan. 2009 glitch/outburst, dur- 
ing the 11 days of segment 3, luminous pulsed 
hard X-ray emission was detected up to ~ 150 
keV. In the next time segment (#4) the maximal 
flux was reached at 70 ± 30 days after the glitch 
epoch, and a gradual decrease by more than a fac- 
tor 10 followed over a period of ~ 300 days till the 
pulsed hard X-ray emission became undetectable. 
This spectral evolution follows the appearance and 
gradual fading of the hard X-ray pulse centered at 
phase ~ 0.59 in the pulse profile (Sects. \4~A\ I5.il 
and FigsM GH>- 

5) During the periods with significant pulsed 
hard X-ray emissions, the total pulsed spectra over 
the broad energy range 2.5 - 150 keV can be sat- 
isfactorily described by a combination of a soft 
power-law (dominating below 10 keV) plus a cut- 
off power-law model (describing the hard X-ray 
emission). The latter component of the spectral 
fits reaches maxima in luminosity (in uF u ) for en- 
ergies in the range 90 - 105 keV and the power-law 
photon index varies with energy from ~ at 30 

keV, 1.2 at 70 keV to 2 around 90 keV 

(Sect.\5J$. 

The actual extrapolation of the spectral shape 
to higher energies (sharp break, bend?) is not con- 
strained by our data. 

6.6. Evolution of the total spectra for X- 
ray energies 1 - 300 keV 

1) Making use of the imaging capabilities of 
INTEGRAL ISGRI, JEM-X and Swift XRT, we 
derived the spectral evolution of the total emis- 
sion of IE 1547.0-5408 over the broad energy range 
1-300 keV starting 2 days after the Jan. 2009 
glitch/outburst till April 2010 (Sects. UTM HOI and 
Fig.EB). 

2) The total spectra can be satisfactorily de- 
scribed by a black-body plus a single power-law 
model over the full 1-300 keV band. The black- 
body component was found to have kT values 
varying in the range 0.57-0.74 keV. The photon in- 
dex of the power-law component exhibited a hard- 
ening with time and a clear correlation with flux 
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in the 20-150 keV band, namely, its value gradu- 
ally changed from ~ — 1.4 directly after the glitch 

in Jan. 2009 to 0.9 in Febr.-April 2010, while 

the 20-150 keV flux decreased by a factor ~ 3 
(Sect.\5^ Tableland Fig.[Tty. 

3) The total soft X-ray flux (1-10 keV) mea- 
sured with Swift XRT (contributions from the 
black-body plus power-law components) is max- 
imal during the first INTEGRAL observation two 
days after the glitch and decays similarly to what 
we measured for the count rates below 10 keV with 
RXTE PCA (Sect.\5J\ FigsEM and\T^. 

4) The total hard X-ray flux (20-150 keV) mea- 
sured with INTEGRAL ISGRI is also directly 
maximal in the first observation 2 days after the 
glitch in Jan. 2009, decaying in flux by a factor ~ 
3 till Febr.-April 2010, when the source was still 
detected in the sky map at a 9.5<r level in the 20 
150 keV band (5ect.E3 Figs. UM CG3 and\T$. This 
behaviour is in contrast to what we measured for 
the evolution of the pulsed hard X-ray component: 
the latter reached maximal flux 70 ± 30 days after 
the glitch, subsequently decaying to undetectable 
levels in ~ 300 days, a decrease in flux by more 
than a factor of 10 (Figs.\Q and \ll\i 

7. Discussion 

The above summarized results on the evolution 
of the high-energy characteristics of IE 1547.0- 
5408 over a period of 27 months, adressing the 
timing parameters, pulse morphologies, total and 
pulsed spectra in the soft and hard X-ray band, 
provide important constraints on the theoretical 
modelling of the production scenario's in the mag- 
netosphere. Particularly interesting is that for the 
first time from the start of an outburst, accom- 
panied by a timing glitch, not only the evolution 
at soft X-ray energies could be studied in detail, 
but also at (non-thermal) hard X-ray energies, and 
the behaviours of distinctly different components 
compared. The measured timing glitch in Jan- 
uary 2009 is itself already interesting, so far being 
the most extreme instantaneous frequency deriva- 
tive jump Ai>/z> of —0.69 ± 0.07 detected, with 
Az> = +(1.30 ± 0.14) • 10" 11 Hz/s. 

We scrutinized earlier publications on magne- 
tar (SGR/AXP) timing, mainly based on RXTE 
observations, and focussed on frequency derivative 
jumps Ai>/v with negative signs of order unity. Al- 



though large grad ual v fluctuations are rep orted 
for SGR 1900 +14 (IWoods et al.lll999l l2002h . and 
SGR 1806-20 (IWoods et al.ll2002l l2007h ~no con- 
vincing sudden instantaneous v jumps are de- 
tected. A possibly similar event as observed 
for IE 15 47.0-5408 may hav e occured for SGR 
1900+14 (|Woods et alJ l2002h near its Aug. 27, 
1998 giant flare of size Az>/;> ~ —0.28, but due to 
a ~ 80 days data gap in the RXTE timing observa- 
tions prior to the giant flare a gradual accelerated 
spin-down could not be excluded. A radiative 
(pulsed and persistent emission) outburst, accom- 
pa nied with a timin g glitch, has been reported 
bv lKaspi et all (l2003h for AXP IE 2259+586, but 
the measured Ai>/z> jump of +1.11 has positive 
sign, and the absolute value of the Ai> jump is or- 
ders of magnitude smaller than that observed for 
IE 1547.0-5408. Also, the glitches with a measur- 
able A;> jump de tected from AXPs 1RXS J1708-40 
(|Dib et al.l l2 008h. IE 1841-045 (iDib et al.ll2008h 



4U 0142+614 (IDib et al.ll2007t 



and IE 1048.1-5937 ( Dib et al 



Gav riil et alj |2011) 



2009) do not show 



the magnitude, sign and relative strength as ob- 
served for IE 1547.0-5408, making the latter glitch 
unique. 

The dramatic decrease in spin-down rate, ob- 
served for IE 1547.0-5408 during the Jan. 22, 
2009 glitch, would for purely dipole spin-down 
naively imply a decrease in effective surface dipole 
field strength by ~ 50%. In the twisted mag- 
netosphere model such a drastic change in v 
and B might be due to a varying twist angle 
with corresponding spectral an d flux changes 
(|Beloborodov k Thompson! l2007j ). The epoch of 
the timing glitch is consistent with the start of 
the outburst, therefore, we assume that the glitch 
triggered the outburst on MJD 54853.035 (TDB), 
Jan. 22, 2009. It is then intruiging that there are 
no signs of a radiative outburst near MJD 55121, 
Oct. 17, 2009, for which our analysis suggests a 
similarly strong Ais/i> jump. 

Recently published works on IE 1547.0-5408 
of which the analyses were performed in paral- 
lel to our work, focussed primarily on the evo- 
lution of the (pulsed) emission below 10 keV in 
time intervals up to ~ 20 days aft er the outbursts 
in Oct. 2008 and Jan. 2009 dNg et all 12011 



Bernardini et al 2011: : IScholz k Kasoi 11201 ill as 



we summarized in the introduction. 

In our work, we revealed in addition the cre- 
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ation of a new transient non-thermal pulse in the 
pulse phase distribution after the Jan. 2009 glitch, 
and could study the simultaneous evolution of the 
pulsed and total high-energy emissions, the lat- 
ter being dominated by an unpulsed component 
(~ 65%). It is striking, that the spectrum of the 
pulsed emission from 2 keV up to 150 keV after the 
glitch is very similar to the shapes of the broad- 
band spectra of the persistent pulsed emissions of 
three AXP s, IE 1841-04 5 4U 0142+61 and 1RXS 



J 1708- 40 (jKuiper et all 120061: Iden Hartog et al 
l2008al |bh. The non-thermal pulsed component 
(20-150 keV) of IE 1547.0-5408 reaches in the time 
interval 15-101 days after the glitch (RXTE time 
segment-4; see Table a flux of (7.5±?;g) • 10~ u 
erg/cm 2 s (see Tabled and an (isotropic) lumi- 
nosity of ~ 13.6 • 10 34 erg/s adopting the d istance 
estimate of 3.9 kpc of Tiengo et al.l ( 2010l ). which 
is about 1.4 times the spin-down luminosity during 
that time segment. There are clearly two maxima 
in the luminosity distribution of the pulsed spec- 
trum, one around 100 keV and one near 1 keV 
(the maximum, likely due to a BB component, is 
reached below our energy interval; earlier works 
did not address the spectra of the pulsed emis- 
sion of IE 1547.0-5408 below 10 keV, nor above 10 
keV). We also see a clear phase separation of ~ 
0.3 between the soft-spectrum pulse and the hard 
non-thermal p ulse, similar to the case of 1RXS 
J1708-40 (e.g. Iden Hartog et al.l l2008bh . In fact, 
also the broad-band total spectrum is reminiscent 
of those measured for these three persistent AXPs. 
However, for IE 1547.0-5408 we are dealing with 
a transient / variable phenomenon above (and 
below) 10 keV, while it was found that the non- 
thermal persistent emission of AXPs appeared 
stable within the statistical errors of ~ 20% over 
more than 10 years. 

The steadiness of t he persistent non-therma l 
emission from AXPs led lBeloborodovl (|2009ll201lh 
to reconsider his scenario's proposed earlier to ex- 
plain the production of this non-therm al compo- 
nent (jBeloborodov &: Thompson! 1200711 . As men- 
tioned in the introduction, he discusses how a star- 
quake can cause convective motions in the crust 
which twist the external magnetic field anchored 
to the surface. The magnetic twist energy is dis- 
sipated over time in the form of radiation, but 
the untwisting occurs in a peculiar way leading to 
the creation of a bundle of electric currents, the 



so-called j-bundle, which shrinks toward the mag- 
netic dipole axis. The high-energy component, 
created by upscattering of thermal X-rays from the 
neutron star surface by an inner relativistic out- 
flow in the quasi-steady j-bundle, is beamed along 
the magnetic dipole axis. The latter can explain 
the relatively narrow high-energy pulses. Further- 
more, his radiative-transfer simulations produce 
spectral shapes that get close to those observed, 
but with possibly too high break energies near 1 
MeV. 

It is very interesting to note that a radia- 
tive outburst accompanied by a strong glitch in 
the case of IE 1547.0-5408 indeed led to a sce- 
nario/geometry which mimics in many respects 
the observational characteristics of the persis- 
tent emission seen from AXPs. But our results 
on IE 1547.0-5408 pose several new observational 
constraints on this model scenario, as well as on 
other theoretical attempts (none of the other mod- 
els addresses the transient hard X-ray (> 20 keV) 
phenomenon of the persistent emission). A con- 
straint on the time scale of creation of the quasi- 
steady narrow j-bundle, is the delay with respect 
to the glitch epoch of 70 ± 30 days, with which the 
high-energy pulse with non-thermal spectrum was 
created. Another important finding is that the to- 
tal non-thermal high-energy emission was present 
immediatly in the first INTEGRAL observation 
starting 2 days after the glitch. This suggests 
that there should be a "corona" around the neu- 
tron star in which non-thermal unpulsed emission 
is produced without delay already a t the start of 
the outburst . The model scenario of iBeloborodov 
(2009|, 2011 ) was proposed to explain the steady 
non-thermal emission from AXPs. What we find 
is a transient pulsed component which decays by 
a factor > 10, becoming undetectable ~ 300 days 
after the glitch. However, the non-thermal to- 
tal (pulsed+unpulsed) emission decayed only by 
a factor ~ 3 over more than a year, seemingly 
stable over the last ~ 100 days. A possibility 
is that a steady state similar to that of persis- 
tent emission of the above mentioned AXPs was 
reached after about one year, but that the pulsed 
non-thermal emission has become too weak to be 
detectable with currently operational high-energy 
instruments. The situation might be representa- 
tive for all those AXPs for which no non-thermal 
emission has been detected (yet), but all have 
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X-ray fluxes in the transition region around 10 
keV that also are weaker than detected for the 
AXPs with reported steady non-thermal persis- 
tent emission above 20 keV up to ~ 150 keV. In 
this respect, we like to note that IE 1547.0-5408 
was in the field of view of INTEGRAL during rev- 
olutions 470-572 (Aug. 21, 2006 - June 22, 2007) 
for an effective exposure of 227.2 ks during a pe- 
riod in which the soft X-ray flux state increased 
at least 16 x its histor ical flux minimum re ached 
in July- August 2006. iHalpern et ail (|2008h con- 



cluded that the source was then recovering from an 
X-ray outburst in between August 2006 and June 
2007. However, we did not find evidence for an 
outburst of non-thermal emission above 20 keV 
from IE 1547.0-5408 in the INTEGRAL ISGRI 
sky maps. The la upper limits to the flux mea- 
surements are at the level of the positive flux mea- 
surements in Revs, 899-912 (around March 2010) 
in Fig. [T31 but still significantly higher than the 
low flux upper limits from the very deep summed 
exposure before 2006 (Revs. 46-411 in Fig.H3]l. 
Also in this case we cannot excl u de th e scenario 
proposed by iBeloborodovl (|20Q9l . 2011 ) , because 
the total and particularly pulsed emission from 
AXPs and also from SGRs is (too) weak for easy 
detection by presently operational instruments. 

We noted the importance of studying the char- 
acteristics of magnetars over a broad energy in- 
terval. In the case of IE 1547.0-5408 we con- 
cluded that discusssions based on spectral model 
fits to the narrow en ergy range 1-10 keV (e.g. 
Scholz fc Kaspi 1 120111 ) leads to apparent charac- 
teristics (e.g. power-law component with spec- 
tral softening with decreasing flux) which have 
to be revised when the broad energy range is 
considered, including the hard X-rays (power-law 
component with spectral hardening with decreas- 
ing flux). Furthermore, broad-band studies of 
the magnetar class are now becoming possi ble, as 
has been published bv lEnoto et al.l (|2010bl ) using 
Suzaku observations of a sample of eight magne- 
tars, including IE 1547.0-5408 in outburst. In that 
paper interesting correlations have been studied 
for the total emission of the hardness ratio, defined 
as the ratio of the flux of the hard, non-thermal 
component (power- law fit; 1-60 keV) and the flux 
of the soft component (a black body, a comp- 
tonized black body or a combination of the two; 1- 
60 keV) with derived characteristics as e.g. char- 



acteristic age and magnetic field strength. Such 
a study ignores, however, that the origins of the 
actual unpulsed and pulsed emissions appear to 
be different (location in the magnctosphcre, pro- 
duction scenario) as is evidenced by the different 
spectra of the pulsed and total (pulsed+unpulsed) 
emissions (particularly above 10 keV), the differ- 
ent pulse profiles (shape and phase) for soft and 
hard pulsed components, and the different evolu- 
tions of all contributing components. Theoretical 
modelling of the physics taking place under the ex- 
treme conditions in magnetospheres of magnetars 
is extremely complex. It is encouraging that the 
amount and detail of observational constraints is 
increasing significantly over the last years. 
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Fig. 1. — The evolution of the total flux (top panel) and pulsed flux (bottom panel) of f E f 547.0-5408 
from Oct. 3, 2008 till December 25, 20f0 as measured by the PCA/PCU-2, summing the signals from all 
three detection layers for the pulse-height-channel range 4-27 (^2—10 keV). We have subtracted the f997 
pre-outburst reference level (see section 13.11 ) for the total rate (top panel; dashed line) and superposed the 
best fit decay model as a solid line (top panel). The middle panel shows the frequency measurements (data 
points), three incoherent fit segments (solid plus dashed lines (= la uncertainty)) and five bold segments 
representing reliable coherent timing solutions. Time intervals for which detailed PCA pulse morphology 
and spectral studies were performed are labeled at the bottom of the middle panel (11 segments). At the 
top of the middle panel 7 intervals are shown for which INTEGRAL observations are available (see Table 
• The dashed line in the bottom panel represents the best fit of a linear decay model for the period MJD 
54868-55555. 
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Fig. 8. — The evolution of the morphology of the (PCA) X-ray pulse profiles in 3 different energy bands 
(see headings at the top of the figure). Time-segment identifiers are shown on the right-hand side, and are 
defined in Table [SJ In particular, segments 1 and 2 refer to the time period between the Oct. 3, 2008 and 
Jan. 22, 2009 outbursts. Segment-3 represents a time period of only 11 days duration starting just after the 
Jan. 22, 2009 outburst, during which rapid and drastic morphology and flux changes occur. 
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Fig. 12.— INTEGRAL ISGRI significance map for the 20-150 keV band of a 20° x 10° field centered on 
IE 1547.0-5408 based on observations performed during INTEGRAL revolutions 899-912 (Feb. 23 - Apr. 3, 
2010), more than a year after the Jan. 2009 outburst. The effective exposure on IE 1547.0-5408 is ~ 717.7 ks, 
and its detection significance is still 9.5c. The location of another radio-loud magnetar, PSR J1622-4950, not 
detected by INTEGRAL, is also indicated by a yellow circle. Other high-energy sources detected by ISGRI 
during revolutions 899-912 are indicated either by name or number. Those labelled with a number are: [1] 
IGR J16207-5129, [2] AX J1619.4-4945, [3] IGR J16318-4848, [4] AX J1631. 9-4752 and [5] 4U 1636-536. 
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